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I. SIMIARY 

To investigate steady-state combustion and to evaluate energy 
sources for sustaining combustion instability, shock wave propagation 
through the combustion zone of a liquid propellant rocket motor is 
under investigation in a theoretical-experimental program. Six equa- 
tions with nine unknowns can be written for the region of interest 
within the chamber. Evaluation of parameters immediately before and 
immediately after the passage of the shock front will be utilized in 
the equation solution. The goal of the research is to produce a method 
by which most of the steady-state combustion parameters can be calculated 
after having measured only the static pressure immediately before and 
immediately after the passage of the shock and the shock velocity. (III-A) 
The experimental program is being conducted with the square-motor 
hardware. An upgrading of the experimental technique from that previously 
reported has been necessary in order to meet the accuracy requirements 
for parameter measurement. Different transducers have been chosen, to- 
gether with vastly improved recording capabilities. The initial shock 
pulse has been tailored to more closely resemble a step-function to aid 

in the determination of the peak pressure. Static pressure survey tech- 
2 

niques have also been further refined. Streak photography is providing 
important secondary data as well as supplying a means for checking the 
prime source of shock velocity measurement, the axial array of pressure 
transducers . (III-B) 

Mixture ratio surveys of the injector types employed in these 
square-motor tests have been conducted with thermocouple probes. This 
technique, using heated pseudo-oxidyzer and room temperature RP-1, has 
been found very helpful in following the extremely rapid mixture ratio 
gradients immediately downstream of the injector face. Spray interpene- 
tration can be modeled via Gaussian distributions. (III-C) 

1 2 

Extensions to the work on droplet population density variations, ’ 
as a possible source of intermediate frequency instability, have continued. 
The effects of additional parameters have been investigated and a non- 
dimensional empirical relationship was derived . (III-D) 

Studies within a pseudo rocket environment to investigate the vapor 


* Refers to Chapter and Section. 
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displacement mechanism for combustion instability have also continued. 

The droplet generator when run in conjunction with the horns, which 
generate the tangential wave motion, has allowed sufficient separation 
of the individual droplets to view the vapor sheaths individually. 

However, the data are qualitative and indications are that significant 
vaporization increases are to be found at locations where velocity or 
pressure is maximized, with the former location producing the greatest 
effect. (III-E) 

A simplified nonlinear model of the gasdynamic oscillations in 
acoustic liner cavities and orifices, and near lined surfaces has been 
developed and analyzed. The resulting nonlinear system of ordinary 
differential equations has been numerically integrated for certain com- 
bustion chamber oscillations, thereby determining orifice velocity and 
cavity pressure versus time. This information is used to calculate the 
liner admittance coefficient which is the important stability parameter. 

(IV- A) 

To evaluate the analytical predictions of the behavior of acoustic 
resonator cavities, a chamber section employing 53 individually adjust- 
able resonators has been fabricated and tested on the 9-inch diameter, 
transverse mode hardware. The tests have involved unlike doublets and 
like-on-like doublets of low and high impingement angle as well as the 
use of an annular motor in addition to the tests on conventional hardware 
where pulsing can be supplied. In these tests there are injector config- 
urations which were "triggered** into high frequency instability and 
such instability could be suppressed with what were often relatively 
small increases in damping. Damp times were influenced by the niimber 
of "active" resonator cavities. The placement of the damping was also 
shown to play a significant role. On the other hand, suppression of 
pressure oscillations in spontaneously unstable tests using unlike injector 
designs would appear to be a more difficult task. Amplitude decreases 
were found but stability could not be attained, at least within the limita- 
tions of the acoustic liner design tested. In those tests the tuning of 
the cavities switched the instability to another tangential mode rather 
than providing stability . (IV-B) 

Extensions of the sensitive time-lag theory to the nonlinear regime 
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have been developed in both the longitudinal and transverse cases. In 
the longitudinal case, both concentrated and distributed combustion are 
considered. Oscillations with shock waves are found, although in a 
limited range of the ^ , n plane near the linear stability limits con- 
tinuous wave forms are found. The possibility of triggering is indicated 
for operation near the stability limits. Distributed combustion is more 
stable than concentrated combustion. (V-A) 

The transverse mode was studied in the special case of thin annular 
chambers. Again, oscillations with shock waves are typical and trigger- 
ing is possible. However, the stability of this mode is not as sensitive 
to the combustion distribution as the longitudinal mode is. In another 
analysis of this transverse case, a simplified droplet burning model was 
employed instead of the sensitive time-lag model. Due to the simplifi- 
cations, no phase lag in the combustion response appeared and only one 
combustion stability parameter (analogous to n ) appeared. (V-B) 

A more complete study of the linear combustion response with a 
quasi-steady droplet burning model has been performed. A small perturba- 
tion analysis has been performed and both the steady-state and perturba- 
tion equations have been solved by numerical means. From these results 
the combustion response factor is calculated. Under special conditions, 
the response factor maximizes in a certain frequency range indicating a 
time-lag behavior. (V-C) 

The transverse wave motion and shock formation in chambers are being 
studied by numerical means. In the first attempt, thin annular chambers 
have been considered. A new numerical technique which involves a rotating 
frame of reference has been employed and the results compare favorably 
with those obtained by a standard technique. In the thin annular case, 
wave steepening and shock formation occur. (V-D) 
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II, INTRODUCTION 

This is the Seventh Yearly Report on the research at Princeton 
University on the problem of combustion instability in liquid propel- 
lant rocket motors. These investigations are being sponsored by the 
National Aeronautics and Space Administration under NASA Contract 
NASr-217 and, in the past, under NASA Grant NsG-99-60. The nonlinear 
aspects of the problem have been emphasized in both the theoretical 
and experimental portions of the program. Presently, the theoretical 
research is concentrating on studies of possible mechanisms of combus- 
tion instability such as droplet burning, the nonlinear gasdynamic 
oscillations associated with combustion instability involving shock-type 
waves in both the tangential and longitudinal modes, the effect of 
damping devices on stability with particular emphasis on nonlinear 
effects associated with acoustic liners, and fundamental studies of the 
basic energy sources. In connection with the theoretical program, an 
experimental program of equally broad range is being conducted. This 
program includes basic experiments with rocket motors and '*cold*' cham- 
bers which have been designed to study various phases of the combustion 
process and/or the instability phenomena to provide understanding which 
would prove useful in the construction of theoretical models. It also 
includes studies of the effects of exhaust nozzles and damping devices 
upon the instability as well as controlled testing of liquid rocket 
motors with parametric variations intended to determine the factors which 
strongly influence the incidence of nonlinear combustion instability. 

In the theoretical study of mechanisms of combustion instability, 
the greatest difficulty involves the construction of the model or models 
which contain the essential physical aspects of the phenomena. In this 
respect the basic experiments are extremely useful. The theoretical study 
of the nonlinear gasdynamic oscillations presents another type of diffi- 
culty. Generally, mathematical techniques for the treatment of nonlinear 
partial differential equations are not well-developed and, in most cases, 
they are undeveloped. Therefore, the theoretician must develop techniques 
as he works with the physical problem. In Sections V-A and V-B of this 
report, certain nonlinear problems are discussed together with the special 
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techniques developed for their solutions. It is the evolution of 
analytical techniques which controls the rate of progress with the 
nonlinear problems. Often, linear mathematical approaches have proven 
useful in effecting a better understanding of the nature of the insta- 
bility. Therefore, a linear treatment has also been used in this report 
to analyze droplet burning. 

As is the norm for any status report, a considerable amount of 
the material discussed herein is in the state of active investigation. 

Thus, only tentative conclusions or observations are sometimes possible. 

At least a mention of one example of each type of work being conducted 
has been made so that this report would prove to be of maximum benefit 
to those readers working in similar areas. 

During the past year the following technical reports, papers and 
book have been (or are about to be) published: 

”Axial Mode Shock Wave Combustion Instability in Liquid Propellant 
Rocket Engines", C, E. Mitchell, Department of Aerospace and 
Mechanical Sciences Technical Report No. 798, July 1967. (NASA CR- 72259) 

"A Preliminary Study of the Effects of Vaporization and Transverse 

Oscillations on Jet Breakup", J. A. Newman, Department of Aerospace 

and Mechanical Sciences Technical Report No. 785, July 1967. (NASA CR- 72258) 

"Experimental Investigations of Liquid Propellant Combustion 
Processes Using Streak Photography", R. M, Williams, Jr., Department 
of Aerospace and Mechanical Sciences Technical Report No. 800, 

July 1967. 

"Acoustic Liner Studies", W. A. Sirignano, L. Crocco and D, T. Harrje, 

ICRPG 3rd Combustion Conference, CPIA Publication No. 138, Vol. 1, 

February 1967, pg. 581. 

"A Theoretical Study of Nonlinear Combustion Instability in Liquid 
Propellant Rocket Motors", B. T. Zinn, Ibid, pg. 444. 

"Axial Energy Distribution Studies in a Liquid Propellant Rocket 
Motor", D. T. Harrje and L. Crocco, Ibid, pg. 523. 

"Longitudinal Shock Wave Combustion Instability in Liquid Propellant 
Rocket Engines", C. E. Mitchell, L. Crocco, W. A. Sirignano, 

Ibid, pg. 395. 

"Effect of the Transverse Velocity Component on the Nonlinear Behavior 
of Short Nozzles", L. Crocco and W. A. Sirignano, AIAA Journal , 

Vol. 4, No. 8, August 1966, pg. 1428. 


8 


Behavior of Supercritical Nozzles Under Three Dimensional 
Oscillatory Conditions , L. Crocco and W. A. Sirignano (in the 
process of being published as an AGARD Monograph). 


Another paper is scheduled to be presented in the next few months 
the title is: 


"Nonlinear Wave Phenomena in Liquid Rocket Combustion Chambers", 

L. Crocco, C. E. Mitchell and W. A. Sirignano, to be presented 
at the Colloquia on Gas Dynamics of Explosions, Brussels, Belgium, 
September 1967. 

Additional background and history of the research on combustion 
instability at Princeton University may be found in the partial list of 
References 1-11. In addition, Reference 12 represents the final report 
in an extensive series of reports on earlier work primarily concerned 
with the linear type of combustion instability (i.e,, spontaneous buildup 
from steady-state operation). 

It should be acknowledged that this work made use of computer 
facilities supported in part by National Science Foundation Grant NSF-GP579. 
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III. ENERGY SOURCES 

A. USE OF SHOCK WAVES TO STUDY COMBUSTION PARAMETERS 

It has long been felt that information on the steady-state combustion 
parameters can be gained by studying the propagation into the combustion 
chamber of an externally generated shock. ^ ^ Consider the following 
sketch: 



If a steady-state is assumed to exist between the injector end (station 0) 
and the station of interest (station 1), the mass, momentum, and energy 
conservation equations between these two sections can be written in alge- 
braic form. The same three principles can also be applied through the 
shock, thus relating the combustion parameters at the station of interest 
just before the arrival of the shock, to the shock velocity and to the 
same parameters immediately after the shock front. Thus six equations 
become available and, as it will presently be seen, they contain the fol- 
lowing nine unknowns: density, pressure and velocity of the gaseous phase, 

and density and velocity of the liquid phase at the station of interest 
just before the arrival of the shock, density, pressure and velocity of 
the gaseous phase at the station of interest immediately after the shock 
front and the shock velocity. Thus if any three of the above nine param- 
eters are measured the other six can be calculated. The most important 
element of this approach is having considered as unknowns the parameters 
immediately after the shock front instead of the parameters way down stream 
after the shock front in, or at the end of, the expansion wave. 

Since the velocity immediately after the shock front is subsonic with 
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respect to the local speed of sound, the shock velocity and all the 
shock parameters are functions of whatever adjustments the gaseous 
and liquid phases undergo in the expansion following the shock front. 

Thus all these adjustments must be included in the equations if the 
parameters in, or at the end of the expansion wave are considered as 
unknowns. This requires that the nature of these adjustments be postu- 
lated and that new unknowns be introduced with the final result that 
too many variables must be measured to solve for the others or that only 
qualitative information can be gathered from the equations. On the other 
hand, the parameters immediately after the shock front can be related to 
those immediately before the shock front by the Rankine-Hugoniot equations 
without knowledge of the down stream adjustments. It is still necessary 
to postulate what changes occur through the front of the shock but these 
changes can be understood more easily than the entire adjustments occur- 
ring in the expansion waves. For instance, considering that the shock 
thickness is of the order of a few mean-free-paths , it can be postulated 
that no droplet shattering will occur at the shock front. The problem 
is that of a shock going around a spherical obstacle. The unsymmetrical 
load will deform the drop which eventually will break, perhaps after some 
oscillations. This is particularly evident if the shock is of moderate 
intensity in which case the unsymmetrical load due to the shock is coun- 
teracted by the inertia and surface tension of the drop. Thus it can be 
postulated that there is no droplet shattering through the front of a 
moderate shock. 

More specifically, it is assumed: 

1. Steady-state between the injector end and the station of 

interest . 

2 . One-dimensional flow. 

3. Instantaneous vaporization of the oxidizer and distributed 
vaporization of the fuel. 

4. That mass, energy, and momentum of the droplets do not 
change through the shock front. 

5. That the shock front is a plane one and moving parallel 
to the flow direction. 

6. That the fuel droplet temperature is equal to the boiling 
temperature of the fuel. 
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7 . That gaseous fuel and oxidizer react immediately to give 
equilibriim reaction products. 

Of the above assumptions only the first one is essential to the appli- 
cation of this method. The others could be partially or totally removed 
by modifying the equations presently used. 

The main goal of this research is to produce a method by which most 
of the steady-state combustion parameters can be calculated after having 
measured only the static pressure immediately before and immediately after 
the shock and the shock velocity. This method of investigation could be 
applied to most rocket engines of current interest. 

The main difficulty, besides that of solving a system of several 
nonlinear equations, is in determining the accuracy of the calculated 
parameters as a function of the accuracy of the measurements and of the 
asstimptions made. 

The Equations 

The equations are now given and then their meaning and their solution 
will be discussed, c.g.s. units are used unless otherwise specified. 
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where are polynomial curve fits of molar specific 

heats at constant pressure versus temperature 
for each of the 8 products 


Equations 25 to 36 are formally equal to Equation 13 through 24 but have 
K , ^ ft / instead of f, ^ B, . 

In the above equations the subscript 0 refers to the injector end, 1 
to the station of interest immediately before the shock arrives, 2 to 
the station of interest immediately after the shock arrives, and ^ to 
the liquid phase (fuel droplets). 

Equation 1 expresses the mass conservation between station 0 and 
station 1 at which some liquid fuel ( ^ = grams of liquid fuel per unit 
volume of the combustion chamber) may exist and move with the velocity 
Equation 2 expresses the momentum conservation between stations 0 
and 1. Since fuel and oxidizer are usually injected at an angle with 
respect to the axis of the combustion chamber {9~) , only the components 
of their velocities along the axis appear in the momentum equation 
( These velocities are related to the feed line velocities 

and pressures by Equations 8 and 9 where a discharge coefficient ) 
takes into account possible losses through the injectors. The initial 
momentum of the liquid, although small, is still important in the momentum 
equation since it is of the same order of the loss of static pressure 
( fo' h ) therefore it is equally important in influencing the veloc- 
ities of the gaseous and liquid phases ( ^ ^^) • 

Equation 3 expresses the energy conservation between stations 0 and 
1. In this equation E; is the internal energy per gram of products (gas- 
eous phase) and H"; is the heat released by the chemical reaction per gram 
of products. In both cases the reference temperature is 298^K . This 
equation was obtained just on physical reasoning but can also be derived 
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from the usual steady one-dimensional conservation equation as follows. 
Starting from: 

^ - XT, 

where (T-29^)^h = latent plus formation enthalpy per unit mass. 

Neglecting viscous dissipation and applying it at stations 0 and 1 it is 
found: 


ojp [ (V - +■ + 

(Toj^ -t- - 2^9) “ 

- -f- A, (f^)^ + r*",: -ff 9/>y] " 


j £ 0 y ^ J ^ 

In the above equation, “/: ^ ^ <p j enthalpies per gram of 

fuel, oxidizer and products respectively. Hence 

X (T, - 29^) = H"/ ^ ^ “ fe //l) enthalpy per gram of 

products. Then comparing this equation with Equation 3 it is seen that they 
are equal if: 

- f,«. A (t, = e,»,A Z /f“ ~ ((if -fe «,A; - H°P 


That the above equation is verified can be seen by substituting in its right 
end side the definitions of Hi and 3- from Equations 10 and 23 and 
using the mass conservation (Equation 1) i.e.: 


- K-.AW, . - )f«,A [- 
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and substituting the following expressions: 


0 ^ 

^ enthalpy of formation of the products per mole of fuel 

(4p ^ enthalpy of formation per mole of fuel 

0"^ o 

= ^<p ^ V enthalpy of formation per mole of oxidizer 
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enthalpy of the products per gram of products 
molecular weight of the equivalent mole 



wT"^ ( a) 


= + "i< “ 


moles of oxidizer per mole of fuel 
grams of products per second 


In the energy equation the term ^A, takes into account heat 

losses through the combustion chamber walls from station 0 to station 1 
(including the injector end surface), and A A (at/^xJ heat losses 
through the plane passing by station 1 and parallel to the injector end 
surface. Finally notice that the heat of vaporization of the fuel, which 
is usually included in the heat of formation of the fuel (since the fuel 
at 298^ K is in the liquid phase) is in Equation 3 considered explicitly. 

Equations 4 and 5 express mass and momentum conservation through 
the shock front. The coefficients / "• ft J ?u takes into account the 
fact that when some liquid fuel is present the effective density of the 
gas is higher. If is the grams of products per unit volume of the 

combustion chamber, the grams of liquid fuel per unit volume of the 

combustion chamber, the grams of liquid per unit volume of liquid 

(equal to the specific gravity of the fuel in the c.g.s. system) then the 
effective density of the products is Pf ^ ^ ^ //L ) • This corrective 
factor is always small ( ^ could possibly be neglected. 

Equation 6 expresses the energy conservation through the shock front. 
The terms E, ^ Ez Hz which are energies and enthalpies per gram of 

products before and after the shock, are related to the temperatures and 
pressures before and after the shock by Equations 10 through 24 and 25 
through 36 respectively. They are also related to the presence of liquid 
fuel by Equation 10. Equations 10 through 24 and 25 through 36 express 
chemical equilibrium among the products before and after the shock. Equa- 
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tion 10 takes into account that not all the injected fuel is available 
for chemical reaction since part of it is still in the liquid phase. 

Thus the fuel to oxidizer ratio changes with the distance from the 
injector end. 

Equation 11 is not an algebraic equation to be solved with the 
others but merely defines which products are being considered. 

Equation 12 defines the molecular weight of the ^Equivalent Mole: 

" which is the mass of reactants per mole of fuel. 

Equation 13 defines the total number of mole per mole of fuel. 

Equations 14, 15, and 16 express the conservation of carbon, hydrogen, 
and oxygen respectively. 

Equations 17 through 21 define the chemical equilibritom between 
the species indicated on the right of the same equations. 

Equation 22 is the perfect gas equation of state. 

Equations 23 and 24 define H“, and C/ in terms of the chemical 
composition of the products and of their temperature. 

As previously stated. Equations 25 through 36 are formally equal to 
Equations 13 through 24 but they refer to quantities immediately after the 
shock. 

All the constants appearing in Equations 17 through 24 either come 
directly from ”JANAF*' tables or are curve fits of data available there. 

Solution of the Equations 

In the first six of the 36 equations presented, the unknowns are: 
f, ^ ' ft, ' A ’ parameters 

are either known constants or are defined in terms of the previous ones 
by Equations 7 through 36. In particular this is the case for the internal 
energies ( , £l) and the enthalpies ( ,//*), which are functions of 

, 7"2 ) ^nd through the chemical processes occurring 

in the combustion chamber described by Equations 7 through 36 . Thus all 
the 36 equations (actually 35, excluding Equation 11) should be solved 
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simultaneously. After a few substitutions they can be reduced to about 
20 nonlinear complicated equations which are still too many to be solved 
simultaneously. Thus they must be broken into groups and iterations 
must be used among groups. In the present attempt to solve those equa- 
bions the fact is used that the product decreases constantly 

with increasing distance from the injector end. It represents the mass 
flux of liquid fuel which must decrease since the liquid fuel can vapor- 
ize but it cannot condense. Moreover this product is equal to zero far 
from the injector end where no liquid is present. Starting where 
Ug A =0 and using experimentally measured and V , 

Equations 10 through 36 can be solved for E, ^Hi ^ ^ for various 

temperatures: X ^ ^ • Curve fitting these solutions the following functions 
are determined: E, = E, ( Tj ; H,» H, (x) ; )^ = X (t;J ; ^ B, ( T^J ; 

Hi = = X, ( ) which give C; ^ H, for any values of 

T, and Ti . The use of these curve fittings is consistent with the fact 
that the chemical equilibrium equations (Equations 17 through 21 and 29 
through 33) and the equations for c^ij. are curve fits of experimental 
results. All these curve fits are expected to be accurate to the fourth 
significant digit. Next two specific volumes are introduced: 


y, ( S.3t^3J ,o) T, 
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Then Equations 1 and 2 are solved for and 2 (^ , and Equations 4 and 

5 for U-/ and while Equations 3 and 6 are left unchanged: 
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When Ft.~ ^C~ ^ Equations 37 and 38 give ‘'7 

and ^ 2 . > Equation 41 gives 'U-, so that Equations 43 and 44 become two 

equations for and to be solved simultaneously. Then Equations 

39 and 40 must confirm that and are null. If they are differ- 
ent from zero, their values are used to calculate ^ then to 

calculate new E", = ^i(t.) ; K. (T.) ) V,-X(^i) / ^i.~ j 

and to repeat the process. There is no guar- 
antee that this iteration on will converge. The fact that the 

behavior of A is known should help. If it doesn't, the problem 

can always be solved for several values of fy and the value which 

satisfies the equations best can be considered the solving one. Preliminary 
results are encouraging but they are not sufficient to deserve discussion. 
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NOMENCLATURE 

atoms of carbons in the fuel molecule 

2 

combustion chamber cross sectional area (cm ) 

2 

chamber walls area up to the station of interest (cm ) 

atoms of hydrogen in the fuel molecule 

specific heat of the liquid fuel (erg/g^K ) 

" " oxidizer " 

” ” " ” gaseous fuel 

” ” ** ” " oxidizer ** 

" " " ** individual combustion products (erg/mole ^k) 

atoms of oxygen in the fuel molecule 
combustion chamber diameter (cm) 

internal energy per gram of products (ref. 298 *k) before arrival 
of the shock (erg/g) 

internal energy per gram of products (ref. 298^^) after arrival 
of the shock (erg/g) 

wall function coefficient 

heat per gram of products (ref. 298 "k) released before arrival of the 
shock (erg/g) 

heat per gram of products (ref. 298^k) released after arrival of the 
shock (erg/g) 

heat transfer coefficient to the walls (erg/sec cm^) 

2 

injector end static pressure (dyn/cm ) 

2 

static pressure at station of interest before arrival of the shock (dyn/cm ) 
M II I! I! I! II after ” ** ” ** 

2 

fuel feeding line pressure (dyn/cm ) 
oxidizer 

heat of vaporization of the fuel (erg/g) 
heat of vaporization of the oxidizer (erg/g) 

distance of the station of interest from the injector end (cm) 

temperature at the station of interest before arrival of the shock ( ®k) 

H II II II II II after ** ** ** ** ’* 


vaporization temperature of the fuel (*K) 

** " oxidizer (^ic) 

"Cp = feeding line fuel temperature (®*f) 


To^ = ” " oxidizer temperature (®*^) 
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V/ = 


v = 

IT = 

<r, = 
= 

Uip = 
= 

WTr = 
WT^ _ 
WTff,= 

X, = 

y, = 

y.= 

^ = 

y® = 
= 

A = 
f, = 

^ = 

n = 

/I = 


velocity of the liquid fuel at the station of interest (cm/sec) 

’’ ” *’ gaseous product at the station of interest before 

arrival of the shock (cm/sec) 

velocity of the gaseous product at the station of interest after 
arrival of the shock (cm/sec) 

feeding line fuel velocity (cm/sec) 

” " oxidizer velocity (cm/sec) 

injector end fuel velocity (cm/sec) 

injector end oxidizer velocity (cm/sec) 

component of t/op in the axial direction (cm/sec) 

It 11 Uo^ *» " 

shock velocity (cm/sec) 

specific volume of the gaseous products at the station of interest 
before arrival of the shock (cm^/g) 

specific volume of the gaseous products at the station of interest 
after arrival of the shock (cm^/g) 

fuel mass flow (g/sec) 

oxidizer mass flow (g/sec) 

molecular weight of the fuel (g/mole) 

*’ " II 11 oxidizer (g/mole) 

equivalent mole (g/mole) 

number of moles of product ^ per equivalent mole 

total nxmber of moles of product at the station of interest before 
arrival of the shock 

total number of moles of product at the station of interest after 
arrival of the shock 

effective moles of oxidizer per mole of fuel 
discharge coefficient 

angle between the injector axis and the combustion chamber axis 

heat transfer coefficient of the gaseous products (erg/sec cm ^ ) 

mass of liquid per unit volume of combustion chamber at the 
station of interest (g/cm^) 

density of the combustion products at the station of interest 
before arrival of the shock (g/cm^) 

density of the combustion products at the station of interest 
after arrival of the shock (g/cm^) 

specific gravity of the fuel (g/cm^) 

" oxidizer (g/cm^) 
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B. EXPERIMENTAL SHOCK WAVE STUDIES 

Associated with the analytical approach to investigate energy 
sources, as described in the previous section, is an experimental 
program being conducted with the square-motor rocket hardware. In 
recent months this program has centered on experimental refinements 
deemed necessary to provide the required measurement accuracy. These 
efforts cover a range of experimental improvements: the pressure 

transducers to be used in future tests are ablatively protected, flush- 
mounted piezoelectric (quartz) units with natural frequencies in the 
400 KHz range, the tape recording system employs equipment with less 
than + l>u, second time displacement between channels and with flat FM 
frequency response to 80 KHz, the streak film recording techniques 
have been further refined together with a new streak window, the initial 
shock pulse has been modified to more closely resemble a step-function, 
and improved means of accurately determining the static pressure profile 
have been developed. Other tests such as the mixture ratio gradient 
measurements described in Section C of this Chapter have also helped to 
better define the early combustion conditions present in the rocket 
chamber . 

Concerning the choice of transducer, the selection was primarily 
based on prevention of any compromise in the range of frequency response 
of the recording equipment (80 KHz). Therefore, the choice was limited 
to transducers whose natural frequency was 3-5 times this value, or 
greater than 240-400 KHz. 

We note with great interest the path pursued by Nicholls et al 
at Michigan in the devlopment of a piezoelectric transducer using a ferro- 
electric element. In these past months we have also been looking into 
the possibilities of producing similar pressure measuring instruments, 
roughly equivalent to the Kistler 603A unit, which would possess charac- 
teristics that might even better fit our individual needs from the stand- 
point of natural frequency, signal output, freedom from acceleration 
effects, etc. This effort is continuing. 
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Currently, the plan is to rely on the Kistler 603A units using 

an .125 inch thickness of RTU 580 ablative coating (a measurement 

16 

technique devloped by Rogero et al at JPL ) . We have found this 

approach extremely useful in our acoustic liner studies. As in the 

previous tests, which were conducted to survey the varying conditions 

13 

that are associated with the passage of a shock wave, six transducers 
have been suitably placed in the square-motor. 

To maintain time agreement between the six transducers to toler- 
ances of less than one microsecond, all transducer signals pass through 
a single recording head onto one-inch tape. Detailed tests on our recently 
acquired Honeywell 7600 equipment indicate that the dynamic time displace- 
ment between outside tracks (worst case) is less than + one microsecond. 
Static time displacement error between any two tracks is of the same 
order. The compound error from static and dynamic combinations also 
appears to fall within a two microsecond band. 

A signal to noise ratio of 48-51 db has been determined using the 
double-extended range (this is the worst case from a S/N standpoint) , 
recording at 120 ips, playback at 1 7/8 ips, and with phase lock on the 
capstan. With 1 7/8 ips playback and phase lock on the tape the S/N ratio 
is reduced to 46 db, however, the time base is preserved. At all higher* 
playback speeds the S/N ratio is not deteriorated by phase locking on 
the tape. These tests were made on worn Memorex tape to duplicate end 
of tape life conditions and hence represent a limiting condition. 

These data indicate that the noise levels represent less than one 
percent of the recorded signal. The overall rise time from 10 to 907o of 
a square wave signal is less than 5 microseconds. These characteristics 
together with excellent time correlation between channels are deemed 
necessary to make the detailed wave shape and location measurements which 
are the key to this experimental program. 

Section A shows that the six conservation equations involve nine 
unknowns, thus for the system of equations to be determinate three of these 
must be measured experimentally. Two sets of these measurable quantities 
seem to present themselves. These are the measurements of static pressure 
behind and in front of the shock, and the shock velocity. Also possible, 
would be the measurement of gas velocity behind and in front of the shock. 
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and the shock velocity. The first set seems to be the easiest and 
will be done with the new pressure instrumentation. The second, which 
makes use of streak photography, involves difficulty in measuring the 
velocity immediately behind the shock and will only be used as a rough 
check on the first technique. However, the streak film can also be 
used for performance measurements since the gas velocity is directly 
related to the characteristic exhaust velocity. Also qualitative 
results should be obtainable from comparison of the pressure traces and 
the streak films, as well as data on the mixture ratio striations at 
various axial positions. 

The principle modifications to the research motor are the addition 
of an 8- inch window section for streak photography and the use of the 
recently developed pulse gun. The 8" window section has 5" slits (of 
1/4 inch width) in each side and utilizes plexiglas windows. Thus, 
streak photographs can be taken with either the combustion light or with 
backlighting. This chamber section is capable of being placed anywhere 
within the first 40" of the rocket motor. Mounted on top of the window 
section are positions for three pressure transducers so that direct 
comparisons of pressure records and streak films can be made. The streak 
camera is of the continuous loop variety (as described in Reference 17 ) 
with the shutter synchronized with the pulse gun in such a manner that 
part of the streak record contains the steady conditions and the remainder 
shows the pressure oscillations following the pulse. A typical streak 
photograph utilizing the new window is shown as Figure 1 . 

The new pulse gun design is the result of the desire to obtain an 

accurate value of the static pressure immediately behind the shock. The 

previously used technique produced an exponential decay following the 

shock and even the best transducers can not respond sufficiently fast to 

18 

accurately determine the true peak amplitude in such cases. To alleviate 

this problem a pulse gun has been designed which generates a step- function 

shaped pulse. This allows the transducer sufficient time to read the 

true pressure level behind the initial shock wave. This pulse is gener- 

* 

ated in a similar manner as before, however, in the present case a much 
larger charge is used and the combustion products are fed into the chamber 
at a controlled rate lasting several milliseconds. A picture of a typical 

* A rapidly burning charge of gunpowder bursts a calibrated disk supplying 
the shock wave. 
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pulse shape can be seen in Figure 2 . 

A modification to the run procedure is a result of the new pressure 
instrtimentation. This limits the possible run time to the order of a few 
hundred milliseconds so that the transducers are not unduly influenced 
(or damaged) by the heat transfer environment. In order to make runs 
for these short times a timer was designed to operate the research motor 
which includes firing the pulse gun and activating the photographic 
system. Since the start transient on this hardware is very reproducible 
and is less than 100 milliseconds in duration, it is expected that this 
system will actually improve data quality by making the operation more 
uniform to run. 

Finally, since steady-state data involving pressure and velocity 
is important in the model, the previous approach using a ten- tube mercury 
manometer has been further improved. A remotely operated ten-passage 
valve has been designed to trap the mercury columns at individual levels 
once chamber pressure has been established. This eliminates the objec- 
tionable column surges(present even with shubbers and porous filters in 
the lines) associated with starting and shutdown transients. Also photo- i 
graphic and parallax problems are no longer a factor. It is unlikely, 
due to valve actuation times, that these data can be obtained during a 
standard run. However, test pressures and mixture ratios can be duplicated 
within a few percent and these static pressure profiles obtained on longer 
duration tests. 

Preliminary tests are now being performed on the overall system 
following initial delays with the Honeywell 7600 recorder which was par- 
tially destroyed when dropped from a cargo plane at Denver. 
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C. MIXTURE RATIO GRADIENT MEASUREMENTS 


The injector is one of the most critical elements in a rocket 
propulsion system. Injectors must be designed to meet a number of 
objectives, including performance maximization, avoidance of unstable 
conditions, and appropriate heat transfer environment ^Unfortunately , 
the mechanisms which link injector design to these objectives are in 
many cases little understood. 

Although obviously the ideal testing situation would have been 
in an actual rocket engine, the difficulties associated with measuring 
mixture ratios under such conditions prohibited this. One must be con- 
tent to assume that the mixture ratio in cold flow tests approximates 
that of actual hot tests. 

Special distribution studies of mixture ratio are of particular 
importance because mixture ratio distribution has been shown to have a 
major effect upon the performance of a rocket engine. Recent investi- 
gations at Aerojet-General Corporation have employed this effect in 

20 

theoretical calculations of performance of rocket engines. This tech- 
nique has resulted in a considerable improvement over the old method of 
calculating performance where this factor was not taken into account. 

Therefore, a method of predicting accurately the spacial distri- 
bution of mixture ratio would be quite useful. Injectors with low 
performance potential could then be identified in the design phase at 
a great saving in time and money which would have been spent in fabri- 
cation and testing. It should be emphasized that the approach described 
here provides data on sharp mixture ratio gradients which are overlooked 
with more conventional sampling approaches. 

The investigations described in this section were meant to shed light 
on this question of how the propellants mix. A grossly simplified model 
has been applied to help explain the observed data and to provide clues 
toward solving the more complex problem of predicting mixture ratio dis- 
tribution for any given injector. 

As suggested before, unstable conditions which may occur in a rocket 
engine are quite often linked to injector design. The same injectors 
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(and other similar designs) which are the topic of this report have 

and are continuing to be used in hot firings here at Princeton to 

1 13 

investigate the problem of combustion instability. ^ 


Experimental Apparatus and Testing Procedure 


The experimental apparatus can conveniently be divided into two 
parts: 1) the testing facilities and 2) the measuring and recording 

apparatus. Further, two radically different measuring and recording 
systems were used at separate times to measure the mixture ratio under 
the face of the injector. The testing facilities include the testing 
booth, the injector and chamber section, the associated systems which 
feed the injector, and a barrel to collect the propellants as they 
drain off underneath. The apparatus employed in the primary method of 
measuring and recording data consisted of a moving thermocouple probe 
rake operated by a motor, Bristol recording charts equipped with re- 
transmitting slide wires, a helipotter, and the IBM 1620 Computer. 

The second method of collecting data, which was used only to check the 
accuracy of the results of some of the tests using the primary method 
(see p^ge 28 ), employed a collector device consisting of 81 glass 
tubes. We shall now discuss each of these systems in further detail. 


Testing Facilities 


The testing booth is approximately eight feet tall with a base 
about four feet square. The injector and chamber section is mounted in 
the top of the booth, and the measuring equipment under the chamber 
section. One side of the testing booth is glass so that the run may be 
observed. The barrel to collect the propellants is located in the lower 
half of the booth. 

The propellants, RP-1 and ethylene glycol (pseudo-lox) are contained 
in separate tanks located on the roof of the building. Both of the tanks 
are nitrogen pressurized and have a capacity of forty gallons. The tanks 
have a pressure limit of 630 psi. Cavitating venturis have been placed 
in both of the propellant feed lines to aid in maintaining the desired 
mass flow rates. 
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Primary Method of Measuring and Recording 

The primary method of data collection utilizes a moving proble 
rake on which ten copper-constantan thermocouples are mounted at 
half inch intervals. Only six of these thermocouples were used in 
the tests, however, since the other four do not pass under the spray 
from these injectors. An ice bath is used for the thermocouple ref- 
erence junction. 

Prior to making a run, the oxidizer is heated up to 180 - 190^ 
Fahrenheit by a cal rod heater located in the roof tank system which 
constantly circulated the ethylene glycol to equalize the temperature 
throughout the tank. The fuel remains close to 60^F. (depending on 
the ambient temperature, storage time, etc.). Each run consisted of 
three different flows while the thermocouple probe moved under the 
spray. The first flow was fuel only, the second flow oxidizer only, 
and the third flow consisted of both the fuel and oxidizer. Knowing 
the temperature of the fuel alone, oxidizer alone, and the mixture at 
a specific point under the injector face, one can determine the mixture 
ratio based on the specific heats of each fluid. 

In Figures 3 - 4 we can see the path that each of the thermocouples 
follows relative to the injector. These paths are perpendicular to the 
rows of fuel elements and oxidizer elements. (In Figure 3 each fuel 
element is represented by and each oxidizer element by "0".) The 
mixture ratio attained by this method along each of the thermocouple 
paths at each of the three axial distances is also plotted in Figures 3 
and 4 . 

Because of the limited number of recording charts which existed until 
recently, and lack of sufficient cables leading from the test cell to the 
recording room, it was possible to obtain readings from only three of 
the thermocouples at one time. Hence, the thermocouple wires were connec- 
ted to a system of relays and wired such that either the odd thermocouple 
readings (i.e. numbers 1, 3, and 5) or the even thermocouple readings 
(i.e. numbers 2, 4, and 6) would be sent over the lines to the recording 
room depending on whether the relay switch was on or off. 

The system of testing was as follows: for each flow the thermocouple 
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probe rake rests under the right chamber edge (in Figures 3-4 ) until 
the propellant injector temperature reaches a steady value. Then the 
probe moves from right to left (in Figures 3-4 ) during which readings 
from the odd thermocouples are being sent to the recording room. When 
the probe reaches the left inside chamber edge an automatic relay 
switch is activated so that readings from the even thermocouples are 
sent to the recording room instead of the odd, and the probe retraces 
the path. A helipotter is used to record the position of the thermo- 
couple probe rake. 

From the recording room the readings from the thermocouples and 
helipotter are sent via the retransmitting slide wires to the IBM 1620 
Computer. Here the data is punched directly onto computer cards during 
the run at the rate of two data cards/second or four data points/second. 

By sending the data directly to the computer rather than taking the 
readings off the Bristol carts the results have been shown to be more 
accurate than hand reduction as well as taking much less time. 

Second Method of Measuring and Recording 

An alternate method of finding mixture ratio data was to use a 

device which collected samples of the spray underneath the injector. 

The collecting apparatus consisted of 81 glass tubes ^ inch in diameter 

arranged in nine rows of nine. The spacing between adjacent tubes was 

% inch. Since the two propellants used are immiscible they separate 

after being collected, providing individual mass flow measurements. 

From these samples the mixture ratio can now be determined as well as 
21 

local mass flow. 

This method was used to check the accuracy of the results obtained 
from some of the tests which employed the primary method of measuring 
and recording the mixture ratio. To do so, the collector apparatus was 
put at the same distance from the injector face as the thermocouple probe 
and the alignment was such that a row of collecting tubes coincided with 
each of the six thermocouple paths. 

The primary method of finding mixture ratio has a great advantage 
over the second method because we obtain a continuous mixture ratio pro- 
file along the thermocouple paths instead of only averaged values at 
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half inch intervals. Since the mixture ratio may vary a great deal 
in a half inch space, the results obtained from the primary method 
are considered more representative of the actual conditions. On the 
other hand, the second method is basically more accurate since mix- 
ture ratio is measured in a much more direct fashion. The second 
method also provides information regarding mass distribution, which 
cannot be obtained from the primary method. 

The IBM 1620 Computer and Calcomp Plotter were used to reduce 
and plot the data respectively from these runs. 

Discussion of Accuracy 


There are of course many sources of error associated with finding 
the mixture ratio using the thermocouple probe. In this section the 
more important sources of error will be discussed and estimates given 
of their magnitude. 

The largest source of error results because the thermocouples 
do not respond instantaneously to temperature variations. As the thermo- 
couple probe moves under the spray the temperature is continually varying, 
and before a thermocouple can become adjusted to a temperature at one 
point it has moved to another point at a different temperature. The 
concentration of liquid streams and droplets plays an important role in 
the thermocouple response. 

The error introduced due to the delayed thermocouple response varies 
in magnitude depending on the degree of temperature variation, the probe 
speed, and naturally is a function of the thermocouple junction size. In 
this probe .010 inch diameter wires were used to maintain position under 
the loads imposed by the liquid streams and droplets. When the tempera- 
ture across the injector face is approximately constant (as it is when 
the fuel alone is flowed or the oxidizer alone) the error introduced is 
small. Where the mixture ratio varied greatly, however, (for instance 
the runs with the Type IV injectors at the 3-inch axial station) quite 
large errors resulted from delayed response. In the majority of runs the 
error introduced by the time delay was estimated as less than + 15%, but 
in one or two cases of extreme mixture ratio variation the error slightly 
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exceeded 50%. 

As mentioned before, the thermocouple probe is placed so that it 
rests just inside the edge of the spray during the initialization (i.e. 
until injector temperatures become steady) and does not go outside the 
spray during any flow. The technique was employed to avoid the larger 
errors that would have resulted from the delayed response if the probe 
had moved outside the spray to inside it or vice-versa. 

One technique which could have alleviated the error due to delayed 
response would have been to drastically cut down the thermocouple probe 
speed (about % inch/second in our tests). Unfortunately, this was not 
feasible because had this been done there would have been insufficient 
propellants to insure a complete run. (The greater part of the propel- 
lant consumption is involved in initializing.) In the future it will 
be possible to take measurements from all thermocouples at once, elim- 
inating three traverses of the probe rake, and thus allowing slower probe 
speeds to be employed. 

Error also arose in determining the position at which the mixture 
ratio is taken. Errors arising due to the helipotter method of determin- 
ing position are no more than + .05 inches. However, there are also 
errors due to alignment in setting up the thermocouple rake or the collec- 
tor apparatus before each run which may reach + .1 inches. 

One error in mixture ratio determination that might also be expected 
would result from air from outside the spray circulating into the spray 
and influencing the temperature measurements. This problem was encountered 
in individual doublet spray tests without chamber wall confinement. However, 
if this were the case in these tests, mixture ratio readings near the edge 
of the chamber would be expected to show greater disagreement with collector 
readings than those in the center of the chamber. This effect was not 
observed. 

One must also consider the accuracy of the Bristol charts, retrans- 
mitting slide wires, and computer operations. This includes error due to 
interpolation of data points from the three flows (fuel, oxidizer, and 
mixture) since the positions of the data points measured (from the computer 
output) during the three flows do not coincide exactly. Inaccuracy in 
mixture ratio due to this factor was less than + 3%. 
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It is interesting to compare the accuracy attained by having 
the data punched up directly on the IBM 1620 Computer with the accu- 
racy that would have been attained by reading and recording the data 
points off the Bristol charts. One run was reduced by both methods 
in order to obtain a measure of this. Coming directly into the com- 
puter is definitely the more accurate of the two methods. The major 
discrepancy between the two methods turns out to be in determining 
position, rather than mixture ratio, since correlating the position 
readings from the chart with the thermocouple readings (on different 
charts) was an error-prone procedure. Position measurements by the 
two methods agreed within + .1 inch and mixture ratio readings within 
+ 2/4. Coming directly into the computer resulted in mixture ratio/ 
position profiles that had a smoother appearance than by the other method. 

Results and Conclus ions 


The experimental findings of this research are summarized in 
Figures 3-4. Figure 3 shows the mixture ratio variation for the 
Type II, 4x4 injector at the three axial stations, while Figure 4 
show the results at the three axial distances for the Type IV, 4x4 and 
Type IV, 6 X 6 injectors placed side by side. 

Basic Assumptions of Model 

1. The differences in mixture ratio distribution between cold flow 
tests and hot firings with the same injector and propellant flow rates is 
negligible, especially near the injector face. 

2. There is no interaction among the various sprays from the dif- 
ferent elements. 

3. Chamber walls have an effect upon mixture ratio distribution 
by reflecting mass from the walls (important for small engines). 

4. All oxidizer elements (they all have like geometry) have the 
same mass distributions. This applies also to all fuel elements. 

5. The mass distribution resulting from a single row of either fuel 
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elements or oxidizer elements along a line perpendicular to the said 
row and parallel to the injector face is symmetric with respect to 
this row. 

Discussion of Model 

The method used in applying the proposed model is straightforward. 
The first step is to make use of assumptions 2, 4, and 5 to set up mass 
distribution profiles along the thermocouple path due to each fuel and 
oxidizer row separately. Then the mass distributions due to the fuel 
rows are added together, as well as the mass distributions due to the 
oxidizer rows. Finally, to each of these sums the effect of mass reflec- 
ted off the chamber walls (assumption 3) is added in so that the total 
mass distribution along the thermocouple path due to each propellant is 
now known. From this the theoretical mixture ratio along the thermocouple 
path is found. By a judicial choice of the magnitude of the key (a, <T) 
one tries to obtain theoretical results that agree well with the experi- 
mental findings . 

The above description leaves two questions unanswered. They are: 

1) What shape shall be chosen for the mass profiles due to each fuel 
and oxidizer row separately? and 2) How can the mass reflected off the 
chamber walls be accounted for quantitatively in the model? The author 
chose solutions to each of these questions which seemed to be the most 
plausible available. The approximate nature of the technique chosen to 
answer these two questions must be kept in mind. As it turns out, dif- 
ferent methods than those chosen would have worked also. A discussion 
of how both of these questions were answered follows. 

It as assumed that the mass profiles due to each fuel and oxid 
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row separately could be represented by the formula f(x) exp 

where a , ^ , and /a. are constants. In the normalized form, the coef- 
ficient a = 1, and then we have the so-called "Gaussian" or "Normal" dis- 
tribution. The mathematics of the Gaussian curve can be found in any 
textbook on statistics (e.g. Ref. 22 ). 

The Gaussian curve exhibits a characteristic bell shape and is 
symmetrical about the mean as is required by assumption 5. Studies made 
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earlier by the author with the collector on single fuel and oxidizer 
elements very similar to the ones employed in these investigations 
indicated that the mass profile along any line parallel to the injector 
face could be fitted by a Gaussian Distribution quite well. It is 
thought that the Gaussian curve could also be applied approximately to 
rows of fuel and oxidizer elements as has been done in the proposed 
model. It turns out that one can manipulate the shape of a Gaussian 
curve quite easily to fit the needs of the model by changing the stan- 
dard deviation (i.e.CT' ) in the preceding formula. 

The problem of accounting for the mass reflected off the chamber 
walls is complicated because one cannot assume that this mass rebounds 
elastically. The spray hits the wall in the form of tiny droplets which 
will be deformed and may break up upon contact with the wall. Some of 
the fluid may in fact run along the wall and not be reflected at all. 

Thus, accounting for where the mass hitting the wall goes becomes an 
exceedingly complex problem. 

For the purposes of the proposed model it was assumed that the 
"coefficient of restitution" for impact of mass in a direction normal 
to the walls is 0.5. This is equivalent to saying that the mass which 
hits the chamber wall rebounds with a velocity (normal component) equal 
to half its incoming velocity (normal component), A discussion of the 
coefficient of restitution (defined for inelastic collisions) may be 
found in certain physics or dynamics textbooks (e.g. Ref. 23 ). Although 
one would expect the coefficient of restitution to vary greatly (between 
0 and 1) for the reflected mass, it is thought that assuming an average 
value might yield a fair approximation. Differences caused because the 
chamber length does not extend the entire distance to the sampling plane 
are neglected in finding reflected mass. 

It was found that the above model could indeed be applied (in a 
trial and error fashion) to produce mixture ratio variations close to 
those observed. Figures 5 and 6 illustrate how this model has been 
applied to thermocouple number three at the three inch and seven inch 
axial stations for the Type II, 4x4 injector. Figures 7 and 8 
illustrate how this model has been applied at the three inch axial station 
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to the Type IV, 4x4 injector and the Type IV, 6x6 injector respec- 
tively for thermocouple number three. The scale for mass flux in these 
figures was chosen arbitrarily. A slight displacement correction has 
been applied to the experimental findings in Figures 7 and 8 where 
careful analysis of the data indicates that the thermocouple probe rake 
was misaligned to the injector. 

Indeed, the proposed model is flexible enough that mass distribu- 
tion curves other than those which were chosen could have been employed 
to fit the V ta equally well. For instance, other Gaussian curves but 
with different standard deviations «r) and which were multiplied by 
different constant coefficients (a) would have worked within a relatively 
narrow range of values. In this fashion one can also manipulate the 
ratio of the area under the oxidizer mass distributions to the area under 
the fuel mass distributions without destroying agreement with experimental 
findings. It is also evident that the assumptions made regarding the 
shape of the mass profiles and the distribution of the mass reflected 
off the walls, however correct or incorrect they might be, could have 
been replaced by other schemes that would have worked. 

In the next section it will be seen how the proposed model can be 
used to deal with the many factors involved and to explain many of the 
observed characteristics of the three injectors. 

Wall Effects and Injector Differences 

In the first place it is interesting to note how necessary it is 
to include the effect of mass reflected off the chamber walls into our 
calculations of theoretical mixture ratio. The reason is as follows. 
Consider for the moment Figures 5-7 in which the theoretical mixture 
ratio distribution for the two 4x4 element array injectors is found. 
Comparing the theoretical mixture ratios at the positions where the two 
oxidizer rows are located reveals that the theoretical mixture ratio at 
the second oxidizer row (going from left to right in Figures 5-7 ) 

must be no more than twice the mixture ratio at the first oxidizer row 
if wall mass reflection is not taken into account. This is because at 
the positions of the two oxidizer rows the oxidizer theoretical mass flux 
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readings are the same (assumption 4) but the fuel mass flux at the 
first oxidizer row can be no more than twice the fuel mass flux at 
the second oxidizer row (from assumption 4 and summing principles) . 

Similar reasoning places other bounds on the problem if wall 
mass reflection effects are not taken into account. For instance, 
for the 6x6 element array injector the theoretical mixture ratio 
without wall reflection at the position of the third oxidizer row 
(soing from left to right in Figure 8 ) can be no more than twice 
the theoretical mixture ratio at the positions of the first two oxi- 
dizer rows. Rather than linger on this point for too long, it is suf- 
fice to say that without taking into account the effect of mass reflec- 
ted off the chamber walls good agreement between theoretical mixture 
ratio variation and experimental findings would have been impossible 
in the wall region. This argument helps to support assumption 3 which 
was made . 

As the axial distance from each of the three injectors is increased, 
the observed variation in mixture ratio becomes less severe. This 
observation is correctly predicted by the proposed model, because as we 
the axial distance one would expect the bell-shaped mass dis- 
tribution due to each row to become wider and flatter (which corresponds 
to increasing in the formula for the Gaussian Distribution). The 

model predicts that mixture ratio variation levels off as the bell-shaped 
distributions become wider (compare Figures 5 and 6 ) . 

The major differences between the mixture ratio variation of the 
Type II, 4 X 4 injector and the Type IV, 4 x 4 injector can be explained 
lis^^b of the proposed model. Consider Figures 9a and 9b which show 
these two injectors. One would expect the shape of the mass distributions 
along a thermocouple path due to separate fuel and oxidizer rows to be 
wider for the Type II injector, since in this injector the thermocouple 
paths along the major axis of the elliptical mass distributions due to 
each element, while in the Type IV injector the thermocouple paths are 
at an angle to the elliptical mass distributions. The proposed model 
would take this effect into account by predicting less severe mixture 
ratio variation for the Type II, 4x4 injector than for the Type IV, 

^ ^ iiijcctor. This in fact is supported strongly by the experimental 
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findings (compare Figure 3 with Figure 4) 

Another experimental finding is that the mixture ratio variation 
along the thermocouple paths for the Type IV injectors tends to be 
more severe toward the center of the injector, which is especially 
evident at the 3-inch axial station (see Figure 4). This observa- 

tion can be at least partially explained by the proposed model in the 
following way. Referring to Figure 9b, one would expect the bell-shaped 
distribution due to each fuel row and oxidizer row separately to become 
wider near the walls parallel to the thermocouple paths since the mass 
distribution (in the direction of the thermocouple paths) of the elements 
nearest these walls would widen upon reflection. According to the pro- 
posed model the mixture ratio variation should become smoother as this 
happens. For the Type II injector, on the other hand, one would expect 
much less reflection from the walls parallel to the thermocouple paths. 

This would help to explain why the effect of more severe mixture ratio 
variation toward the center is not observed in this injector. 

Suggested Method of Predicting Mixture Ratio 

In this last section it will be shown how the model proposed in 
this paper might be modified and extended in order to solve the more 
complex problem of predicting mixture ratio variation for a given injector. 
Such a method of predicting mixture ratio variation would be based on 
assumptions 1, 2, and 3 listed previously, and would proceed as follows. 

First, a detailed mass survey using a collector device would have 

to be taken (at the proper flow rate) for each fuel and oxidizer element 

of the injector which has a different geometry. If an element is of the 

like-on-unlike variety, we would also need mixture ratio information. 

Note that it would appear to be necessary to take surveys of these single 

elements at only one axial station. This is possible because it can be 

shown (Ref. 24) that the flow rate along a radial line from the point 

2 

of impingement of a spray is given by where 

and G^ are the local mass fluxes on the given line corresponding to 
the distances and respectively. Thus if the mass distribution 

in one plane due to the element is known it can be found for any other 
plane, however, in the confined space available in the rocket motor tested 
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these regions overlap. 

Second, a method of accounting accurately for the distribution 
of the mass flux reflected from the walls would have to be developed. 

This is very essential to the problem for small rocket motors as the 
investigations described in this paper show. 

Lastly, the mass distributions of each element, along with its 
position and orientation would be plugged into a computer program. 

The mass distributions would then be superimposed and wall reflections 
taken into account by the computer program. Theoretical mixture ratio 
distributions could then be calculated. 

It is evident that if this method of predicting mixture ratio is 
to work, many details must be ironed out in the future. The problem is 
obviously not a simple one especially when combustion effects are present 
in an actual firing. 
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D. DROPLET POPULATION DENSITY INVESTIGATION 

A possible source of energy for buzz-type instability, i.e, , 
injection-coupled phenomenon, is the grouping of the injected droplets 
into high and low density clumps. Viewed from the combustion "zone'* 
these variations are akin to mass fluctuations even though the feed 
system might be supplying a constant mass flow. 

The preliminary results of the study of variations in droplet 
population densities in impinging jet sprays were reported in last year's 
progress report.^ The cold flow study of this phenomenon is now finished 
and is more completely outlined in a MSE thesis to be released shortly as 
a technical report. In the past year, additional variables have been 
investigated with the primary emphasis being placed upon the separation 
of the liquid properties of viscosity and surface tension; until recently 
these properties were combined into one dimensionless parameter. In 
addition to the actual parametric studies an attempt was made at improving 
the existing data recording and reduction techniques. 

The experimental investigation of the droplet population density 
variation involved the following parameters: breakup frequency, chamber 

pressure, injection velocity (A p) , injector diameter (proportional to 
drop size), and the fluid properties, density, viscosity and surface tension. 
Also included in the dimensional analysis (described later) is the gas 
viscosity, although this parameter was not varied in this cold flow study. 

As reported last year, the ratio of the number of droplets in the wave max- 
imum to that in the wave minimum was also measured. This was done as a 
check for the existence and magnitude of possible resonance effects. The 
primary test procedure was to fix all of the variables except chamber pres- 
sure, and to then make a series of four runs at various pressures. After 
data reduction of the streak film, a curve of breakup frequency vs. chamber 
pressure was obtained for the various test conditions. The previously 
reported results ^ were obtained with only two fluids, water and 507o 
water-507o methyl alcohol. The resulting data did not allow the fluid prop- 
erties to be separated so further tests were run with a surface tension 
reducing agent added to water. This allowed the surface tension to be 
significantly varied while maintaining viscosity and density constant. Tests 
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were also conducted with 100% methyl alcohol to more closely simulate 
the actual rocket motors. 

From the resulting data an empirical relation was obtained using 
nondimensionalized analysis. This approach was as follows: 

assuming that frequency = / * ^ 
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From the test data one can now obtain the exponents. These in 

fact are: c^ = -.25, c^ = 3.25, c^ = .8, and c^ = 2.9. From the 

resulting final empirical plot (Figure 10) one also gets the empirical 

1 ^ / 

constant (A = 2 x 10 ) . Thus the final cold flow empirical result is: 


- 2X/0'^ 



Also seen in Figure 10 is a curve resulting from the data obtained 
from transverse rocket tests exhibiting intermediate frequency pressure 

10 

oscillations. As can be seen, this curve has a constant of A = 7.0 x 10 
It is believed that this difference in the empirical constants can be 

attributed to the inability to include the gas viscosity directly in the 

* 

analysis, except in conjunction with the liquid viscosity. The large 
increase in gas viscosity (about two orders of magnitude) in the rocket 
motor seems to account for the large change in the empirical constant. 

Thus, we see that this experimental evidence seems to support the hypo- 
thesis that the droplet breakup phenomenon does affect intermediate fre- 
quency combustion instability. 

It should also be noted that no effects were found to indicate any 
enhancement of the natural breakup phenomenon by operating the chamber 
under oscillation conditions. However, it is felt that the low amplitude 
oscillations induced by the siren (peak- to-peak amplitude approximately 
107 o of steady-state pressure) and the lack of sensitivity of the streak 
technique for such measurements tend to make inappropriate any conclusion 
concerning the possibility of resonance effects. 

To further investigate resonance effects and to reduce the effect 
required to obtain and reduce data, various electro-optical techniques 
were investigated. These included the direct use of photocells and photo- 
transistors to ’’look” into the chamber and the indirect use to reduce the 
streak film. However, after several attempts to obtain a workable system 
these techniques were abandoned in favor of the previous approach. This 
was necessary in order to insure that all data was reduced within a reason 
able time, since it appeared that significant sophistication and hence 

* This parameter could not be significantly varied in the cold flow 
experimental apparatus since most gases have similar viscosities at 
room temperature. 
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further development time was required in the case of electro-optical 
techniques before obtaining a working system. 

Returning to the rocket motor tests displaying intermediate fre- 
quency pressure oscillations, several interesting features are evident 
in the pressure traces (Figure 11 ). These include phase relationships 
between transducers separated by 90^ and phase relationship between the 
maximum and minimum envelopes of any one transducer. Such effects are 
not readily explainable since acoustic propagation should be found for 
any pressure variation and hence such phase differences should be assoc- 
iated with much higher frequencies. 

E. VAPOR DISPLACEMENT AND STREAM BREAKUP EXPERIMENTS 

Another energy source within the rocket motor is the coupling of 

the pressure oscillations via associated vapor displacement with the 

resulting changes in the instantaneous burning rate. Referred to as the 

"vapor displacement mechanism" this coupling is important in the case 

of transverse mode instability. Over the past few years experiments 

on pseudo rocket hardware have sought to add quantitative data to the 

1-3 

qualitative evidence of the importance of this mechanism. A recent 

25 

technical report covers one portion of this investigation in detail , 

In that report experiments performed to investigate the breakup of a 
vaporizing liquid Freon jet are described. Behavior of the jet has been 
studied under a variety of steady and oscillatory pressures. 

It has been observed that chamber pressure in the pseudo rocket has 
no significant effect on jet length in the laminar regime so long as the 
fluid is not rapidly evaporating. For the fluid whose vaporization rate 
is substantial, a decrease in chamber pressure and the subsequent increase 
in vaporization, leads to a slight decrease in jet length. This effect 
becomes more pronounced the larger the jet diameter. In the turbulent 
region the effect of pressure depends on the fluid under observation, the 
jet velocity and jet diameter. For a very slowly evaporating fluid, cham- 
ber pressure always accelerates jet breakup. This can be attributed to 
the increased effect of air friction. The length of a vaporizing Freon 
jet, however, sometimes decreases with chamber pressure and under other 
circumstances increases. Generally, the evaporation process becomes more 
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important than air friction in controlling jet length when the jet 
diameter is large and the jet velocity is high. 

When the gas in the chamber is subjected to acoustic oscillations, 
jet breakup is generally enhanced regardless of whether the fluid is 
rapidly evaporating or not. The effect, however, is only measurable in 
the laminar regime for the range of oscillation amplitudes employed. 

At low steady-state chamber pressures, oscillatory pressure (standing 
mode at a pressure antinode) and oscillatory velocity are both about as 
effective in promoting the breakup of a non-evaporating fluid. At higher 
steady pressure levels, pressure oscillations (standing mode at a pressure 
antinode compared to the spinning mode) appear to be more effective than 
velocity oscillations. The behavior of a vaporizing fluid is essentially 
the same as a non-vaporizing fluid except that velocity oscillations be- 
come as effective in promoting breakup as pressure oscillations at high 
chamber pressure. As chamber pressure is lowered and vaporization rates 
increase, velocity oscillations become more effective than pressure 
oscillations . 

Under none of the circumstances encountered in the pseudo rocket 
experiment did the spinning mode enhance jet breakup to a greater extent 
than the standing mode. 

Viewing the results of droplet spray fan mixing in a hot rocket 
environment, tangential displacements resulting from velocity variations 
were observed. A sodium seeded jet was used as a tracer and indicated no 
detectable unidirectional motion within the range of the test conditions. 

Very recent experiments have concentrated on separating the injected 
Freon jet into droplets via vibrations of the hypodermic tube injector. 

It was found that uniform, well-spaced droplet streams could only be formed 
from specially vibrated laminar jets. Thus we were concerned with low vel- 
ocity, low Reynolds number flow. Through extensive experimentation it was 
found that droplet stream stability required relatively high, yet still 
laminar, velocity; large hypodermic diameter; and high vibration frequency. 
This stability is essential in that photographs of an established stream 
in the various chamber modes must show droplets in approximately the same 
position if there are to be definite picture comparisons. Also, streams 
without good stability are very often destroyed by the violent pressure and 
velocity oscillations that are studyed. Since once the stream was established 
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it could not very easily be checked until after the entire run, much 
time was wasted in taking and developing pictures which had already 
been destroyed. 

In the present study the major interest was in obtaining maximum 
droplet separation (so as to be able to distinguish individual vapor 
sheaths in shadowgraph pictures). This requirement made very low jet 
velocity, and low generator frequency necessary. Clearly then there 
was a compromise between separation and stability; and experimental pro- 
cedure here became a matter of sensitive touch, patience, and often a 
good deal of luck. Figure 12 shows such preliminary data. 

With the Strobotac placed on the spark source side of the chamber, 
a laminar jet of about one-half inch is established. A chamber pressure 
of at least 20 psig was necessary when Freon 114 was used to keep this 
liquid from vaporizing instantly. Laminar flow was easily seen as a 
straight-line jet free of perturbations. With its amplifier turned up 
to full volume, the audio generator was turned on and a stable droplet 
was established at about 100 cps and the drops were made visible and 
stationary by adjusting the Strobotac in its low frequency range. From 
this point both the generator and stroboscope frequency were continuously 
dialed down until a stable stream could no longer be maintained. During 
this slow-down process the jet velocity also had to be cut back. Once 
the optimum separation point was reached the generator frequency was 
raised slightly to insure stability throughout the following run. 

Placement of the neoprene feed tube was crucial - as much of the 
tube as possible had to be suspended so as to feed freely and vertically 
into the hypodermic. If the tube hung from the generator its weight created 
a sideways force which disrupted uniform flow sending droplets in more than 
one direction in a spray. 

Many low-frequency streams required the full power of the 40-watt 
home-made amplifier and care was taken not to overheat the hypodermic 
voice coil. 

The tentative findings of this study is that droplets are more rapidly 
vaporized in velocity as opposed to pressure oscillating fields. 
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IV. DAMPING DEVICES 

A. THEORETICAL STUDIES OF FLOW BEHAVIOR WITH ACOUSTIC LINERS 

One of the essential goals of the study of acoustic liners is to 
determine the optimum cavity design for a chamber oscillation of a 
particular frequency or range of frequencies. For small amplitude oscil- 
lations the cavity resonant frequency is determined by the Helmholz 
resonator formula, but for larger amplitudes there may be a significant 
deviation from this value. One of the purposes of this analysis is to 
try to determine the extent of this deviation. 

The main dissipation mechanism under consideration in an acoustic 
liner is the loss of kinetic energy of the gas jets as they leave the 
orifices. This kinetic energy is converted into turbulence of the gas 
in the chamber or the cavity, and hence no longer contributes to the 
energy of the oscillating wave. This mechanism is nonlinear because it 
involves the square of the velocity. The equations describing the process 
are not solvable analytically and one resorts to other techniques. In 
this analysis the equations contained several parameters, and solutions 
for particular values of these parameters were found by numerical integra- 
tion on the 7094 computer. 

The known quantity here is taken to be the chamber pressure as a 
function of time and the desired results are the orifice velocity and 
cavity pressure oscillations. 

Two sets of equations have been developed, making a different assump- 
tion about the flow in the orifice in each case. Based on the second order 
analysis of Sirignano the assumption was made that the velocity is 

constant along the length of the orifice at any instant of time. A second 
case was run for comparison in which the assumption was that the mass flux 
is constant along the orifice at any instant of time. 

The model used to develop the equations is isentropic flow in both 
directions through the orifice but with different mean gas temperature in 

the chamber and in the cavity. Thus we assume the cavity gas maintains 

y 

the relationship ^ = C where C is a constant in time. The chamber 
gas is also assumed to follow this relationship but with a different value 
of the constant. 
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In both cases the variables were non-dimensionalized with respect 
to the mean chamber conditions. The velocity direction convention can 
be seen from Figure 13 . 

It is also assumed that the acceleration of the gas to the velocity 
u at the entrance of the orifice occurs in a very short distance and 
with no lag in time. Thus, the pressure at the orifice entrance may be 
expressed as 




< 



where the choice of subscripts will depend on the flow direction and where 


W, = 

JTcpfT 

w^ may written in identical form. 

In the first case we assume velocity constant along the orifice 
length and a linear density variation, hence at any instant of time, 
u = const and ^ = b^ + b^x where x is axial distance along the 
orifice and b^ and are constants. 

Under these restrictions, the one-dimensional, integrated momentum 
equation, applied to the orifice as a control volume gives us the following 
result for positive flow. 


diL - -a] 


where 


W/ r: 




and the continuity equation applied to the cavity yields, 




Yu ^ 


Similarly for negative flow, we obtain 


where 


W, n 






and 

In the above, cx = ^JN . Also, pressures are non- 

dimens ionalized with respect to p^ , the mean chamber pressure, orifice 
velocity with respect to , the mean speed of sound, and time with 

respect to Jtl^^ . 

It must be noted that as a consequence of the assumptions these 
equations do not match together at zero velocity, that is 


^ dt 



In the second case, the assumption was that the mass flux is constant 
along the length of the orifice. With this assumption, the momentum 
equation gives, for positive flow 

if- = ^UO-wf)^'K-^] 
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where 




And continuity for the cavity results in, 

^ = Yo&^nn 


While for negative flow. 


if ='> In - K-^] 


where 
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and 




4 - Yd. 

Yck ^ ^ 


The values of w^^ .and must be determined from the following 

equations. 


c/-w‘r-'w, = 






and 




m 




In the above equations, in is nondimens ionalized by • These 

last equations must be solved using an iteration scheme for particular 
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values of m and or . 

In the second case the equations are not discontinuous at zero 
velocity. 

We have neglected such effects as flow separation at the entrance 
to the orifice, nonzero velocity in the chamber and temperature changes 
in the cavity. Further study will have to determine the actual magnitude 
of these effects, but it is hoped that this analysis will be a useful 
first approximation to such a complex problem. 

The scheme chosen for integrating the equations is a form of the 
Milne variable- interval technique. This technique was chosen because it 
allows a smaller interval to be used where it is needed, without taking 
as much computer time as a Runge-Kutta scheme for that same interval size. 
Also there was a well tested program available at the Guggenheim Computing 
Center which utilized this technique. 

Even with the equations in nondimensionalized form there are many 
parameters which must be chosen before a solution can be obtained. These 
involve the particular form of the input (chamber) pressure, the geometry 
of the cavity-orifice combination, and the ratio of the isentropic con- 
stants in the cavity and the chamber. 

The programs are presently arranged to take an input pressure express- 
ed as a Fourier series, but it would be relatively easy to substitute some 
other functional form. Also, the programs were run only for a value of 
the ratio of specific heats (V) of 1-2 but it would be a simple matter 
to change to some other value. 

It was decided that the most convenient form of output would be in 
the form of Fourier components, both since it is compact (yet definitive) 
as well as the fact that one of the ultimate uses of these programs will 
be to calculate values of admittance coefficients for an acoustic cavity. 
The admittance coefficient is determined by the ratio of the complex vel- 
ocity at the entrance to the complex chamber pressure, where the real parts 
of the velocity and pressure are the first harmonic components expressed 
as a cosine with a phase angle. 

A check was made to assure that there was no mean mass flow in either 
direction (this would be inconsistent with our analysis) and it was found 
the mean flow in each case was less than one-half of one percent of the 
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first harmonic mass flow component, which was considered within the 
range of numerical accuracy. 

In order to have a constant reference for the frequency, it has 
been expressed as a ratio of the actual nondimens ional frequency to 
the Helmholtz resonant frequency. We find that. 

Here, f is nondimensionalized with respect to . So in the 

data a frequency of 1.0 corresponds to the Helmholtz resonant frequency. 

All of the values of orifice velocity and pressure amplitude are 
first harmonic values. The chamber pressure wave consisted of a constant 
(of value 1 because of the nondimensionalization) plus a first harmonic 
term only. As this work progresses further, it is intended to investi- 
gate the effect of higher harmonic terms in order to more closely approx- 
imate actual chamber waves. 

The range of interest of the parameter ^ is quite large, for 
we see that the resonant frequency of the cavity-orifice combination from 
the Helmholtz formula is. 


where is the cavity sound speed. Since the resonant frequency desired 
for a wide range of rocket chambers varies over about two orders of magni- 
tude and possible orifice lengths vary over about one order of magnitude, 
the value of ^ could vary over as much as six orders. Due to reasonable 
limitations on computer time we were able to study only two values of 
These two values were 0.1 and 0.01. It can be seen from Figures 14 to 21 
at the lower value of ^ the pressure amplitudes are much lower. This 
seems reasonable in view of the fact that for a given orifice length and 
orifice area, decreasing ^ amounts to increasing the size of the cavity 
which would naturally result in lower cavity pressure amplitudes. Also, 
for the lower value of , the velocity amplitude becomes less sensitive 
to frequency variations, at least up to 1.2 times the Helmholtz resonant 
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frequency. This would seem to follow from the fact that the velocity 
derivative is proportional to the pressure difference across the orifice 
and as the cavity amplitude decreases this difference is less affected 
by phase differences. Hence ultimately the velocity amplitude varies 
less . 

The parameter is determined mainly by the ratio of the tem- 

peratures in the cavity and the chamber. With a little manipulation we 
can see that, taking mean conditions in the cavity and the chamber, 



where the subscript * refers to a mean cavity condition, and stands 

for the molecular weight of the gas. Since we have no idea of what the 

value of Pq/Pv- should be, we cannot pick both it and the value of 

T,/T^ . We are fortunate then that for V = 1.2 this ratio occurs to 
* 0 

the 1/5 power in o< , So, assuming it to be close to 1, we calculate 
from the following, 



As a check on this assumption it was found from the integration 
that the value of Pq/p.,. was generally in the range 0.95 to 1.05, which 
indicates that we have an error on the order of 1%. 

The programs were run for three values of the parameter oC . The 

value oc = 1.0 gives the case in which both chamber and cavity are at 

the same temperature. In this situation the Helmholtz formula should be 

most accurate and this is indeed seen (for low amplitudes) in Figure 14 . 
The value c< = 0.15 corresponds approximately to the case of a *'hot” 
cavity into which hot chamber gases have been pumped. Finally, the value 
c< = 0.05 simulates the ’’cold” condition in the cavity, which is normally 
the situation when instability begins. 

It appears that as the value of o< is reduced from 1.0 the pressure 
peak moves down in frequency and also the velocity reaches higher values 
at low frequencies of the order of 0.2 times the resonant frequency. 
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The nonlinearity of the problem under consideration is very 
evident from the data. We could reasonably expect if the problem were 
linear that we would have maximum cavity pressure amplitude, maximum 
velocity amplitude, and 90 degree phase lag of the cavity pressure 
related to the chamber pressure at the point corresponding to the 
natural frequency of oscillation of the orifice-cavity combination. 
However, the situation in this case is not quite so simple. While 
the velocity amplitude seems to peak right near the Helmholtz natural 
frequency, the pressure peak is strongly affected by the value of the 
parameters oC and ^ . In addition, the cavity pressure peak is 
shifted by an increase in amplitude of the chamber pressure wave. 

As far as the two different cases are concerned, they show the 
same general trends but the amplitudes of the velocity and pressure 
responses differ in most cases. For example, in the velocity response 
for Case 1 the velocity amplitudes around the resonant frequency de- 
crease with decreasing c< , whereas in Case 2 the amplitudes remain 
very nearly constant. 

Further points will have to be taken to determine if there is an 
optimum value of or if the trends shown with the two taken in this 

analysis continue throughout the range of interest. It may also be 
important to determine which case is giving valid results, if either, 
with regard to the velocity amplitude change with decreasing oc . 

Since increasing corresponds to increasing the cavity temperature, 

which would occur if hot chamber gases were being pumped into the cavity, 
an increase in velocity amplitude could mean an increase in damping 
capability in this situation. 
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NOMENCLATURE 

A - orifice area 

Jl - orifice length 

V - cavity volume 

a - local sound speed 

p - pressure 

^ « density 

u - orifice velocity 

m - orifice mass flux/unit area 

C - isentropic pressure density constant = C) 

^ - molecular weight 

Subscript I - condition in chamber 

Subscript II - condition in cavity 

Subscript 1 - condition at chamber end of orifice 

Subscript 2 - condition at cavity end of orifice 

Subscript 0 - mean chamber condition 

Subscript * - mean cavity condition 
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B . EXPERIMENTAL LINER STUDIES 

To supply further data for the liner models and to simulate the 
actual firing conditions experienced by an acoustic liner, an experi- 
mental program has been initiated. The background data for these 
liner studies is the history of both linear and nonlinear high-frequency 

transverse mode combustion instability that has been documented over the 

1-13 

past years at Princeton. The rocket hardware is the 9-inch diameter, 

heat-sink chamber using either 12 or 16 spuds of like-on-like or unlike- 
type doublets arranged on a 7 or 8 inch injection diameter with thrust 
levels ranging from 1000 to 1333 lbs. Chamber pressures have been nom- 
inally set at 150 psia through the choice of nozzle throat size and flow 
rates. Pulsing is achieved through the use of tangentially oriented 
pulse guns. The propellants are LOX/Ethanol. 

Key to the acoustic liner studies is the array of individually 
adjustable acoustic resonating cavities which are arranged circumferen- 
tially in three rows (axially spaced .75, 2.0 and 3.25 inches from the 
injector face) on a four-inch long section of the 9-inch diameter test 
hardware. The test configuration is shown in Figure 22 . Based on the 
environment present within the cavities, the screw- type length adjustment 
is made to achieve the desired frequency of this Helmholtz resonator. 

Each of the individual cavities has six .1 inch diameter orifices, \ inch 
long. The acoustic cavity itself is nominally 3/4 inch diameter and is 
capable of length changes from zero to two inches. The zero length con- 
dition makes the cavity inactive and hence it does not affect events in 
the chamber. This point has been checked (i.e., with all cavities at the 
zero length setting) and no noticeable change in the triggered” first 
tangential mode oscillation amplitudes was observed. 

Of considerable importance in the early tests, as well as in some 
preliminary tests with a single active cavity, were the temperature mea- 
surements which were made to better determine the cavity environment that 
fixes the cavity lengths required for proper resonant tuning. Following 
the temperature variations with .005 inch diameter W 5% Re and W 26% Re 
thermocouples, one finds that, even with stable operation, temperatures of 
1300^F are reached momentarily following ignition and then the temperature 
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exponentially returns to a few hundred degrees above ambient. Ini- 
tially, an unstable run behaves much the same way, however, rather 
than the exponential fall off the pumping action of the hot gases 
causes the temperatures to reach temperatures as high as 2300^F, 

Little variation was found between the conditions in the first row 
cavities versus third row cavities (.75 vs. 3.25 inches from the 
injector face) . 

Tests were initiated on the like-on-like doublet injector con- 
figuration consisting of 16 spuds with .059 inch orifice diameters 

3 

spaced evenly on a 7-inch injection diameter. Previous tests indi- 
cated that this injector, chamber combination was nonlinearly unstable 
to tangential pulses. Table I illustrates this fact, noting that the 
oscillations initiated in the chamber by the pulse gun fell in the range 
of 226-236 psi pk-pk, while the resulting instability maintained a 
200-250 psi pk-pk amplitude. The frequency was in the 2900-3000 Hz 
range and damping did not occur. 

In the first tests, only three resonating cavities were activated. 
These were grouped in a 72^ sector and were located in the first row 
(i.e., .75 inch from the injector face). The noticeable change with 
these cavities present was the approximately 40% reduction in instability 
amplitude . 

Moving the three active cavities to a 120^ circumferential spacing 
(i.e., related to the tip positions of the highly effective 3-bladed 
baffle designs - protection to the 3T mode) caused the resulting IT mode 
instability to momentarily exist at 74-76 psi pk-pk amplitude and then 
damp in 170-285 seconds to combustion noise (c.n.). As in previous exper- 
iences with lower amplitude oscillations the frequency was depressed to 
2720 Hz, a result of reduced gas mixing and a lower speed of sound near 
the injector. 

Increasing the number of cavities at the 120^ location to 3 at each 
point (in three different rows) the damp time was reduced to a narrow 
time band of 28-30 milliseconds. See Figure 23 for damping trend. 

Further proof that the 3 0 120^ (.75”) resonator cavity grouping 
was a marginal damper, was shown by tests at the 2 inch location. Here 
the amplitude stayed at the 78-97 psi pk-pk amplitude level and damping 


never occurred. 
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Reductions of damping times to the minimum value would have 
required tests with all 53 cavities tuned for the 2800 Hz range, 
i.e., the .375 inch cavity length. However, these tests were missed 
in the concurrent survey with the unlike doublets (the tests that will 
be described next) . The prime reason for the concurrent program was 
the relative ease of changing injectors as contrasted to the time 
required to reset and pressure check the cavities when considerable 
numbers were used. However, data on a mixed group of cavity settings 
is shown in Table I. The .141 inch length was designed to suppress 
the 2T mode. Shorter damp times 12-22 millisec. were indicated in 
these tests as was expected. 

The next series of tests were conducted on spontaneously unstable 
injector designs. Here the injector consisted of 12 unlike doublets 
with .082 inch orifices. The propellant choice was again LOX/Ethanol 
primarily to eliminate carbon logging of the acoustic cavity orifices 
as well as to minimize collection of fuel in the cavities and thus 
cause possible difficulties in subsequent firings. 

The stability behavior of the unlike doublet injector design was 
that of spontaneous IT spinning mode instability on every test. The 
final amplitude reached, as noted in Table II, was nominally 340 psi 
pk-pk. The frequency was 3150 to 3185 Hz, an increase over the most 
unstable like-on-like tests. 

Because this injector design was spontaneously unstable and to 
minimize the danger to the water-cooled Dynisco transducers, due to the 
associated high heat transfer environment, the liner studies were initia- 
ted with nine cavities open at the 120° spacing which had proved to be 
dynamically stable for the like-on-like injector. In these tests the re- 
sultant IT spinning mode amplitude averaged 275 psi pk-pk (a 20% decrease). 
* o 

Seventeen active cavities @ 120 produced similar results. Twenty- 

* 

six cavities produced noticeable amplitude reduction with several regions 
of distinctly different amplitudes being observed, a situation possibly 
attributed to the mixed mode behavior. With 52 active cavities tuned to 
the IT mode frequency the mode within the chamber switched to the second 

* Occasionally there were gaps in the orderly arrangement of the 

resonating cavities to allow for chamber pressure transducers and 
the pulse gun. 
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tangential and the amplitude was found to fluctuate between 50 and 150 
psi pk-pk. In an attempt to eliminate both the IT and 2T modes, 27 
cavities were adjusted to resonate at the 2T mode and 26 at the IT mode. 
Again spontaneous instability occurred, with 133 psi pk-pk amplitudes, 
this time returning to the IT mode. 

An attempt was made to tune the cavities for the IT mode based on 
the "cold'* cavity conditions associated with the start of a rocket firing. 
By coincidence this was the same cavity length required for the 2T mode. 
The instability again was spontaneous and amplitudes of 216-258 psi pk-pk 
were recorded. This test serves as an illustration of the effectiveness 
of cavity tuning although it was designed for another purpose. 

Failing to reach dynamic stability with the liner capabilities assoc- 
iated with the 53 cavity configuration, two simple baffle designs were 
explored for instability suppression. Sketched below is a four-bladed 
baffle with hub and a six-bladed baffle. Both were two inches in length 
and fabricated from 1/8 inch thick copper. 



Sketch Showing 
Two Baffle Designs 

The four-bladed baffle design experienced 2T mode oscillations of 
190 psi pk-pk amplitude at 4650 Hz. With only two transducers @ 90 circum- 
ferential spacing it was not possible to differentiate whether this was a 2T 
standing or spinning mode, however, the location just downstream of the 
baffle tips suggest^ a 2T standing with the nodal lines fixed by the baffles. 
The six-bladed baffle limited the instability to primarily the IT mode 
(evidence of higher modes, 2T, 3T, and combined modes IT-lR are found in 
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the frequency analysis) with only 20 psi pk-pk oscillation amplitudes 
(0 only 2460 Hz). 

This later approach was the most successful of those tried to 
suppress the spontaneous -type instability associated with this unlike 
doublet injector design. 

Testing with the acoustic liner was resumed on the annular motor. 

This chamber design retained the 9-inch diameter chamber sections but 
added a 7-inch diameter center body, and replaced the conventional nozzle 
with a 24 multiple orifice design. The cylindrical length was increased 
to six inches (previously a 4- inch long chamber was used). The chamber 
Mach number at the point of nozzle entry was increased from .05 to .15. 
Chamber pressure was retained at the 150 psia level. 

Although the annular chamber provided an opportunity to test the 
acoustic liner at reduced frequencies (1600 vs. 3000 Hz for the IT mode), 
and with shock- type waves, the prime purpose for this hardware design 
was quite different. As outlined in Section B of Chapter V a theo- 
retical model of combustion instability has recently been formulated by 
Professor Crocco which is based on an annular motor configuration and a 
simplified nozzle condition which can be represented by the multiple 
orifice design. The annular chamber will be used primarily in the exper- 
imental checkout gf certain predictions of this theory. 

From a combustion instability standpoint the annular chamber exhib- 
ited a noticeably different behavior than the conventional design. Elements 
that were dynamically stable in the conventional chamber (recessed like-on 
like doublets with 40° impingement angles, again with ,059 inch diameter 
orifices) were found to produce self-triggered instability in the annular 
motor. Nonlinearly unstable elements (the standard like-on-like doublet 
designs) were self-triggered into IT standing mode instability in the 
annular chamber. The spontaneously unstable unlike doublets were not 
tested in this hardware. 

Table III illustrates the stabilizing effect of the acoustic liner 
in the annular chamber with 53 and 26 cavities adjusted to a cavity length 
of 1.15 inches (i.e., the 1600 Hz range). Sources of popping, which are 
believed due to wall splashing and the presence of the centerbody causing 
local propellant accumulations, were present in the damped tests as well. 
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Sxjmmarizing the tests just described - it would appear that for 
this hardware and propellant combination (LOX/ALC) there are configu- 
rations which can be ’’triggered” into high frequency instability and 
such instability can also be suppressed with what are often relatively 
small increases in damping. Damp times were influenced by the number 
of active resonator cavities as shown in Figure 23. The placement of 
the damping was also shown to play a significant role. On the other 
hand, suppression of pressure oscillations in spontaneously unstable 
tests using unlike injector designs would appear to be a more difficult 
task. Amplitude decreases were found but stability could not be attained, 
at least within the limitations of the acoustic liner design tested. In 
those tests the tuning of the cavities switched the instability to another 
tangential mode rather than providing stability. 

The nature of the energy distribution as measured from the radial 
variation in the pk-pk pressure oscillations was also investigated in 
the liner studies. Based on the unlike doublet tests one might conclude 
that the energy loss was primarily in the outer zone. Figure 24 illus- 
trates this point. The instability prior to the addition of the liner 
was found to fit the theoretical distribution^^ at the mid-radius loca- 
tion. Using the same injector, with the acoustic liner in place and all 
the cavities activated, the wall pk-pk pressure measurement was depressed 
507o while the mid-radius point reading was lowered only 32%. 

Complicating the picture, however, were the observations made with 
the like-on-like injector design, which provided wall pressure amplitudes 
close to the case just discussed. These tests also indicated a similar 
non- theoretical profile. Based on this later data it is difficult to 
interpret the results at this time. 

The tests referred to here used a single Kistler 615 transducer in 
the two pressure measurement locations to eliminate any possibility of 
transducer error. Additional run to run monitoring was done with a 
Dynisco transducer located at the wall. Further testing of the radial 
energy profile is anticipated in the near future using the newly acquired 
tape recording system and employing Kistler 603A transducers. 

Another experimental-analytic program, falling within the scope of 
the acoustic liner studies, involves the local phase conditions prevailing 
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between the chamber wall and the acoustic liner cavity. Using the 
model described in the previous section a brief survey was taken as 
to what values of phase might be anticipated under the conditions 
present in a typical rocket chamber. In an H shaped survey as shown 
in Figure 25 several important regions were investigated. First, the 
conditions representing the onset of instability (i.e., a cold resonator 
cavity) were surveyed over a range of amplitudes, to the point where 
orifice choking intervened (M > 6). A similar survey with fully-estab- 
lished instability (i.e., hot cavity) was also made with choking occur- 
ring at a higher value of ( . Linking these two surveys was the trans- 
ient condition of rising cavity temperature. Moving from the cold to 
the hot cavity conditions would predict phase variations from 147° to 
72°, 

Although a number of preliminary rocket tests have already been 
performed to check these predictions it was felt that the data should 
be retaken using the new recording system prior to making final compar- 
isons. Such tests are now in progress. 
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V. EXTENSIONS TO THE THEORY 


A. SHOCK WAVE INSTABILITY IN THE AXIAL MODE 


This topic is the subject of a Ph.D. dissertation that will soon 
be available as Technical Report 798. Consequently the complete analysis 
will not be presented in this yearly report, but rather a summary of the 
principal results of the work will be given. 

In summary, the investigation has succeeded in developing an anal- 
ytical method for predicting whether periodic finite amplitude axial 
oscillations may exist in a given liquid rocket combustion chamber if 
the combustion processes occurring there can be represented using Crocco*s 
time-lag postulate. Both discontinuous (shock type) periodic waveforms 
and continuous periodic waveforms for the oscillations have been found. 

The dependence of the amplitude and the waveform of these periodic oscil- 
lations on n and T , the combustion parameters, has been investigated. 
Two particular combustion distributions have been considered; the first 
being the limiting case of concentrated combustion, the second being the 
case where combustion is distributed throughout the chamber so that the 
steady-state velocity profile satisfies the equation = 0(Kg) . In 
general, it may be said that, in contrast to earlier expectations, ’ 
the possibility of supporting both continuous and discontinuous finite 
amplitude pressure oscillations using a combustion model with a character- 
istic time that is of the order of the period of oscillation has been 
demonstrated. 

Critical to the successful analysis of the problem was the selection 
of an expansion parameter proportional to, and of the same order as, the 
Mach number at the entrance to the supercritical nozzle. If a different 
expansion parameter, based on the difference between energy input and 
withdrawal from the chamber, was employed, it was found that discontinuous 
periodic oscillations could not be predicted when the combustion processes 
were represented by the time-lag model. It is to be noted that this latter 
kind of expansion parameter was used in all previous linear and nonlinear 
work which employed the time-lag postulate. 

The Mach number at the entrance to the nozzle serves as a measure of 
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the total amount of combustion occurring in the chamber in the steady- 
state and also as a measure of the mean flow out of the nozzle under 
steady conditions. It can therefore be thought of as a measure of the 
deviation of the flow field in the chamber from a simple acoustic field 
with no mean flow. Picking the expansion parameter to be proportional 
to, and of the same order as, the Mach number in this analysis therefore 
caused the governing partial differential equations for the flow field 
to take the form of simple acoustic equations to lowest order. No energy 
or mass addition needed to be considered to this lowest order. The com- 
bustion (driving) terms and the nozzle (damping) term appeared only in 
the second order equations and boundary conditions. Because of this, 
standard techniques could be employed in order to relate the first order 
thermodynamic and gasdynamic dependent variables to an arbitrary function 
of order M (M is the Mach number). Application of the solid wall boun- 
dary conditions then showed that this first order function was periodic 
with period equal to the acoustic period for the fundamental mode of oscil- 
lation. A second order analysis then had to be pursued in order the find 
the form that the arbitrary first order function must take. The second 
order nozzle and combustion terms were found to appear only in the boundary 
conditions (concentrated combustion) or in the partial differential equa- 
tions and the boundary condition (distributed combustion) . In either case 
it was these terms, combined with the nonlinear wave propagation terms, 
that determined the form of the arbitrary first order function. 

The second order partial differential equations that included the 
combustion terms and nonlinear effects were solved with relative ease, 
particularly in the case of concentrated combustion. However, the integra- 
tion of the equations introduced two arbitrary second order functions. 

Thus, three arbitrary functions appeared in the second order expressions 
for the dependent variables. Three conditions therefore needed to be 
applied in order to find an equation determining the first order function. 
The three conditions available were the boundary conditions at x = 0 
and X = 1 , and the periodicity requirement. The boundary conditions 
were rather easily applied, though the fact that Riemann invariants are 
continuous to second order in shock strength through shocks of the proper 
family had to be used in order to apply the boundary conditions correctly 
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when discontinuous oscillations were considered. The application of 
the periodic condition, on the other hand, presented some difficulties 
because, in general, the period of the assumed oscillation had to be 
considered different from the period of the first order function. In 
other words, the period could not be taken a priori to be equal to the 
acoustic period but, rather, was taken to have a value within a correc- 
tion of order M of it. The difficulties occurring because of this 
difference in periods were overcome by introducing a stretched time 
coordinate, © , to replace the physical time coordinate, t . This 

insured that to all orders of approximation the period of oscillation 
in the stretched time variable would be equal to the acoustic period. 

When this transformation was made, the periodic condition was implicitly 
applied and, therefore, the application of the boundary conditions to 
the second order solutions in € , x space produced a nonlinear equa- 

tion determining if periodic solutions for the arbitrary first order 
function, f , were possible and, if so, what their form and dependence 
upon n and jx was. ( ^ is the stretched time-lag that enters when 
the time coordinate is stretched.) Also, the equation determined whether 
the periodic solutions were continuous or discontinuous. 

For the case of concentrated combustion the equation determining 
the form of f was a nonlinear, first order ordinary differential equa- 
tion with a retarded variable. For the general case of distributed com- 
bustion the equation was a nonlinear differential-integral equation with 
a retarded variable. However, for the special combustion distribution 
with ^ , the governing equation could be transformed into the same 

equation as the one governing f for the case of concentrated combustion. 

Only a little information on the behavior of the ordinary differential 
equation for the case of concentrated combustion could be found analytically, 
however these findings were of some importance in setting up the numerical 
procedure for the integration of the equation. The most important result 
of the analytical examination of the equation was probably the discovery 
that discontinuous solutions to the equation with the same n value, but 
with jX values displaced symmetrically from the line ^ = 1 , could be 
related by a simple 180^ rotation. Because of the relationship between 
solutions symmetrically displaced from ^ = 1 , curves of amplitude as 
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functions of n and ytt- were found to be S3raunetric with respect to the 
line This meant also that the nonlinear stability limit had 

to be symmetric with respect to the line ^ = 1 . The periods of the 
oscillations that were solutions to the governing equation, on the other 
hand, were found to be antisymmetric with respect to yw- = 1 . in par- 
ticular, solutions along the line yU- = 1 had to be resonant oscillations, 
solutions with jju > 1 below resonant oscillations and solutions with 
^ < 1 above resonant oscillations. These S 5 mimetry characteristics 
of the equation meant that ntomerical integration of the nonlinear ordinary 
differential equation only had to be carried out for ^ ^ 1 in order 
to find all solutions in the interval 0 ^ jll ^ 2 . This interval 
corresponds to the range of interest for the fundamental longitudinal mode 
of oscillation. 


In addition to the symmetry properties of the solutions to the equa- 
tion governing f for concentrated combustion (or distributed combustion 
. , du 

with "5^ = ^ two additional analytical findings are of interest. The 
first of these is the requirement that the integral of f over one period 
of oscillation is zero. This means that the deviation of the mean pressure 
from the steady-state value must be zero to first order, even when discon- 
tinuous pressure waves are present. Moreover, for the case of distributed 

combustion with 4^ = 
dx 

entropy variation is equal to zero, since this quantity is proportional 
to the mean pressure. It should be noted that the comments just made with 
respect to the mean pressure and first order entropy need not necessarily 
be true for general combustion distributions with = O(M^) . However, 
it is likely that, because of the spreading of the combustion, the deviation 
in the entropy from its steady value should be small, and, since the two 
quantities are proportional, the mean pressure change should also be small. 

The second additional analytical result was the finding the discon- 
tinuous solutions existed right up to the stability limit when = 1 . 

8 

This agrees with a result of the analysis performed by Sirignano using 
a combustion model that did not consider phasing effects. He also found 
discontinuous solutions right up to the linear stability limit. Since 
the line A = 1 represents the line along which energy is added in phase 
with the pressure oscillations in the present analysis, the agreement 


, it implies that the value of the first order 
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between the two investigations is obvious. Moreover, it was found 
that, as one moved along the line ^ = 1 , away from the linear sta- 
bility limit, the amplitude of the shocks increased linearly with the 
distance from the linear limit. The same result was found by Sirignano 
using his different combustion model. 

As far as the actual numerical integration of the ordinary nonlinear 
differential equation resulting from the study of concentrated combustion 
is concerned, little can be said without going into great detail. However, 
a method of functional iteration had to be used because of the appearance 
of retarded variable, ^ , in the equation. This method converged rapidly 
to a proper periodic solution, if one existed at all, for the n and 
values considered. This numerical technique was used to find both continu- 
ous and discontinuous solutions. However, a convenient method of finding 
approximate continuous solutions valid close to the linear stability limit 
was also developed. In order to find these approximate solutions the gov- 
erning ordinary differential equation was treated by a Fourier analysis 
technique in which the quantity f/M (or taken to be a small 

quantity that vanished at the stability limit, in agreement with earlier 
linear analyses. The agreement between the approximate solutions and the 
numerical solutions for continuous oscillations close enough to the linear 
stability limit was found to be very good. However, at larger distances 
from the stability limit, considerable divergence between the two methods 
of solution was observed. 

Numerical integration of the equation in question for a wide range 
of values of n and yU- permitted the construction of a nonlinear stability 
limit for discontinuous solutions as shown in Figure 26 . Inside this 
curve it is possible to find stable, periodic discontinuous solutions. 
Outside the curve no periodic discontinuous solutions at all could be found. 
Comparison of the nonlinear stability limit with the linear stability 
limit shows that three regions with different stability characteristics 
exist. The largest region is the region where the area inside the linear 
stability limit is also inside the nonlinear stability limit. Since, by 
linear analysis, the growth of small perturbations was predicted here, 
while, at the same time, discontinuous periodic pressure waves are also 
predicted by the present nonlinear analysis, the conclusion is reached that 
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the predicted discontinuous oscillations in this region are the final 
form that pressure oscillations would take in a rocket engine that 
was intrinsically unstable because of its n and yU, values being in 
this region of linear instability. In other words, the discontinuous 
periodic oscillations are taken to be the regime form of instability 
in this region. It should be noted that much of this region is far 
from the linear stability limit, in regions where f/M can be quite 
large. However, the analysis is ^till as3nnptotically valid for small 
anough M , even when this is the case. In contrast to this, the valid- 

g 

ity of the nonlinear work of Sirignano, which predicted only contin- 
uous finite amplitude oscillations using the time-lag model, was re- 
stricted to regions very close to the linear stability limit. 

The second region of interest was found to be the region that 
is inside the linear stability limit but outside the nonlinear stability 
limit for discontinuous oscillations. The situation here is found to be 
really quite similar to the one just mentioned, except that the final 
form that intrinsic instability takes in this region is found to be char- 
acterized by periodic continuous oscillations of finite amplitude, rather 
than by periodic discontinuous oscillations. 

The final region of interest is the area that is outside the linear 
stability limit (in a region of linear stability) but inside the nonlinear 
stability limit. It can be argued that this is a region where the possi- 
bility of triggering discontinuous periodic pressure oscillations exists. 
This is because, for a given n and yU in this region, two periodic 
solutions are found, (see Figure 27 ), a large amplitude discontinuous 
solution and a smaller amplitude solution that was either discontinuous 
or continuous, depending on the normal displacement of the n , yn point 
under consideration from the linear stability limit. The smaller ampli- 
tude solutions are interpreted as being unstable solutions and thus as 
representing a triggering limit, in the sense that a disturbance with an 
amplitude smaller than that of the small amplitude solution at a given 
n and A would decay to zero, in agreement with the predictions of 
linear analysis, while a disturbance with an amplitude greater than that 
of the small amplitude solution would grow in time and have the large ampl 
tude discontinuous solution as its limit cycle. In other words, though 
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combustors with n and in this region should be intrinsically 
stable, it should be possible to trigger a discontinuous periodic 
oscillation by the introduction of a disturbance with large enough 
pressure amplitude. Note that the final form that triggered insta- 
bility takes is predicted in the present analysis. Earlier work by 

g 

Sirignano predicted that triggering was possible, but, because of 
the restriction of his analysis to continuous waveforms, the final 
form and amplitude of the oscillations so triggered could not be pre- 
dicted . 

In the last two regions of the n , yic plane mentioned, contin- 
uous periodic oscillations were found. In the region where no trigger- 
ing was possible, but where linear instability was predicted, the con- 
tinuous oscillations were found to be stable. In the other region, 
where triggering was possible, the oscillations were taken to be unstable. 

In both cases the question of how well the continuous solutions matched 
the discontinuous solutions arises. In the former situation the matching 
occurs at the boundary of the region, in the latter it occurs inside the 
region under investigation. The correspondence between the amplitude 
and period of the continuous and discontinuous periodic oscillations was 
found to be very good in both cases. The matching of the amplitudes of 
oscillation between the continuous and discontinuous solutions is shown 
in Figures 27 and 28 for two regions on the n , ^ plane where, respec- 
tively, triggering is possible and impossible. The matching of the first 
order periods of oscillation for discontinuous and continuous oscillations 
in the same two regions is shown in Figures 27 and 30 . In both cases 
the matching point occurs at a normal displacement from the linear stability 
limit where the amplitude of g becomes asymptotically equal to Tf . 

This corresponds to the point where the nonlinear ordinary differential 
equation passes from a region where singular behavior is possible to one 
in which no such behavior is possible. The former is true for amplitudes 
of g greater than T* ^ , the latter for amplitudes less than *T ^ 

Therefore, it is found that a natural transition from discontinuous to 
continuous waveforms is inherent in the governing ordinary differential 
equation. 

The wave shapes themselves were found to be strongly dependent on 


67 - 


the parameters n and . This is true for both continuous and dis- 
continuous oscillations. The continuous oscillations exhibited pressure 
waveforms, generally of a sinusoid nature, with a steep, sharp positive 
peak followed by a shallower negative peak. The sharpness and steepness 
of the positive peak was observed to increase (as well as the amplitude 
of the oscillation) as the normal displacement from the linear stability 
limit was increased. The wave had an arbitrarily steep slope arbitrarily 
close to the value of the normal displacement where the amplitude of g 
became equal to ^ ^ . This behavior is shown in Figure 31 . 

The salient feature of all the discontinuous waveforms that were 
found was the negative infinite slope at Jc time units after the shock. 
This can be seen for the particular waveform shown in Figure 32 . This 
reflects the fact that a discontinuity in pressure occurring at a given 
instant has an effect on the mass release and therefore on the waveform 
exactly Xx. time units later. Because the combustion term, or mass 

' o 

generation term, is one order of magnitude smaller (0(M ) vs. 0(M)) than 

the shock itself, the delayed effect is not felt as a discontinuity in f 

but rather as an infinitely rapid continuous change in f . 

Another characteristic of the discontinuous waveforms was their 

change in shape as the normal displacement from the linear stability limit 

was varied. In general, as the displacement from the linear stability 

limit was increased, the waveforms tended toward the classical "sawtooth" 

2 7 

form found in one dimensional cavities with no combustion. As the dis- 

tance from the linear stability limit became small, on the contrary, the 
wave shapes tended to have a more sinusoid nature (except, of course, for 

the infinite slope at U, time units after the shock) , similar to those 

26 ' 

predicted by Chester for driven oscillations near resonance. This 
changing character of the waveforms is shown in Figure 33 . 

In conclusion, it should be emphasized that the possibility of pre- 
dicting discontinuous oscillations which represent the form that regime 
oscillations in a liquid rocket motor with the combustion processes repre- 
sented in a gross way by the time-lag model, has been demonstrated. 
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B. NONLINEAR TRANSVERSAL INSTABILITY 

In addition to the nonlinear study of axial mode instability described 
in Section A, work is proceeding on the investigation of nonlinear insta- 
bility for the transversal modes of oscillation. 

The general mathematical techniques used in the study of longitudinal 
oscillations are also employed in the transversal mode investigation. The 
Mach number is used as the expansion parameter in the power series repre- 
sentations of the dependent variables necessary for the solution of the 
governing system of nonlinear partial differential equations. The combustion 
chamber is taken to be a thin annulus and the oscillations assumed to be 
purely transversal. Under these conditions the lowest order partial differ- 
ential equations are simply the linear equations of acoustics, as was the 
case for the longitudinal mode of oscillation. This equation is solved in 
terms of an arbitrary function representing the fundamental spinning mode 
of oscillation. An equation determining the wave shape of this first order 
function is found by solving the second order partial differential equations 
and by applying the injector and nozzle boundary conditions as well as the 
conditions at the travelling shock wave and the periodicity requirement. 

In order to avoid secular terms in the equations and in order to properly 
apply the periodic condition it is necessary to stretch the independent vari- 
ables in the form of a power series in the Mach number. The expansion of 
the independent variables is analogous to the expansion of the time coordinate 
in the longitudinal study . 

Two combustion models are considered in the case of transversal oscillations. 
The first uses the heuristic representation of combustion provided by the 
time-lag model. The second employs a simplified droplet evaporation model. 

When the time-lag model is used the equation controlling the waveform of 
the first order function for the fundamental spinning mode is of the same form 
as the corresponding equation controlling the form of the waveform for the 
fundamental longitudinal mode. In fact, through the use of a simple transfor- 
mation the two equations can be made identical. Thus all the results found 
for the case of concentrated combustion in the longitudinal mode can be used 
directly for the fundamental spinning mode. These results will, of course, 
include the prediction of the waveforms for both discontinuous and continuous 
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oscillations. It should be remarked that, though a distributed combustion 
zone is considered in the transversal investigations, the fact that the 
addition of mass is distributed in the axial direction and not in the 
transversal direction means that no integral terms appear in the governing 
equation for the spinning mode of oscillation as they did in the governing 
equation for longitudinal oscillations when distributed combustion was 
considered. 

The simplified droplet evaporation model mentioned earlier as the 
second combustion model which is considered is defined by the following 
assumptions: 

1) The transformation into burned gases immediately follows the evaporation 
of the propellants, so that the rate of combustion coincides with the rate 
of evaporation. 

2) The rate of evaporation is determined in a quasi-steady fashion by the 
well known empirical e^ 2 ^ression used by other authors (i.e., Priem and 
Heidmann, NASA TR-R67) with the conduction term neglected with respect to 
the convection term. 

3) The drag of the droplet is considered negligible. 

With a few additional assumptions, the terms can be ordered in such a way 
that for vanishing injection flux the combustion stays distributed and that 
the first nonlinear terms of the gasdynamical equations, as well as the 
effects of the pressure oscillations and of the transversal velocity oscil- 
lations on the combustion rate are all contained in the second order equations. 

The solution then can proceed in a way similar to that used for the 
time-lag model, with the result that instability is possible only if the 
quantity 

^ 
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is smaller than a certain value, under which waves of a very sharp type are 
possible independently of the frequency of the system. 

The fact that the result is independent of the frequency can be attri- 
buted to certain of the assumptions used in the treatment. More sophisticated 
assumptions should bring in evidence the effects of the frequency. 
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C. DROPLET BURNING UNDER OSCILLATORY CONDITIONS 

A linear analysis for the dynamic behavior of droplet vaporization 
process in a rocket combustion chamber with transverse acoustic waves is 
developed. Small perturbation of the oscillation components of all the 
physical quantities is assumed in this study. The dynamic response shows 
a strong dependence on the oscillation frequency. The characteristic 
curve of the response reaches a peak positive value at a particular fre- 
quency and becomes negative when the frequency is larger than a certain 
value. The effects of combustor conditions (pressure, temperature, and 
gas exit velocity) and initial conditions of the liquid drop (temperature, 
size and velocity) are shown. Comparison with the nonlinear numerical 
result is discussed . 

Inside a rocket combustion chamber, acoustic mode instability takes 
place when the energy of the combustion is generated in some particular 
manner that amplifies the acoustic oscillations. A combustion process 
can be considered to consist of several processes in parallel or in se- 
quence: atomization, vaporization, mixing and chemical reactions. Certain 

assumptions and approximations must be made in order to simplify the problem 
for analytical treatment. Droplet vaporization is one of the processes 
to be suggested as rate-controlling. That means, the vaporization of liquid 

drop is the controlling factor for the rates of heat release and mass addi- 

29 

tion. In the work of Priem and Guentert, variations in vaporization 
rate were shown to be able to affect combustion stability significantly. 

30 

This vaporization process was studied more thoroughly by Heidmann and Wieber 
from a nonlinear numerical analysis by using a model predicting the vapor- 
ization rate for the entire drop life. In their work, acoustic oscillations 
of pressure and velocity are superimposed on the mean vaporization process. 
Each drop’s history is considered for different boundary conditions and 
summation of all the drops is made to obtain the time variation of vapor- 
ization rate in the combustion chamber. This nonlinear numerical analysis 
motivates the investigation of a linear approach based on the same combustion 
model. A more rigorous mathematical treatment with the assumption of small 
perturbations is developed for the case with transverse oscillations. 

It is assumed that the mixing and chemical reaction rates are much 
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higher than the vaporization rate. So the combustion rate is vaporiza- 
tion controlled. Several asstimptions are made in the derivation: 

1) The mean steady flow is one dimensional. 

2) The mean combustion gas pressure and temperature are essentially 
constants . 

3) All droplets have the same initial radius, temperature and 
velocity (i.e., at the injection face). 

4) Mass transfer results from a driving force of concentration 
gradient only, 

5) The droplet temperature is considered to be uniform over the 
whole droplet. 

28 

Details of the derivation can be found in the work by Priem and Heidmann. 
The vaporization rate of a single droplet is 


w Mum 

where 


p'' 2.7n. r* D’' 



( 1 ) 

( 2 ) 

(3) 


and Num is the Nusselt number for mass transfer and is related to the 
Reynolds number by the semi-empirical relation: 


M 


am 


2. -V 0.6 Sc^ R 


♦ "k 


(4) 


and Reynolds number is defined by: 


R< 


2 r* p’* \ %S\ 


(5) 



The heat which is transferred from the gas to the droplet is: 

Qt = Z (T‘-Te ■) "Z 


where Nuh is the Nusselt number for heat transfer: 


Nu. 


I \ 


2. + 0. 6 Pr^ R 


is the correction factor for heat transfer. It represents the ratio 
of heat that is conducted with mass transfer to that without mass transfer. 
The value of is always less than unity since the mass transfer reduces 

the heat transferred to the surface of the droplet. 




- lo 


and 


w*G 


Q. 


^ is the heat received by the droplet and causes the droplet heat-up, 

Q: Q* - W* /V* 


The mass of the droplet is: 




7c r 






The vaporization rate per unit volume is: 




n* vj* 


The equation of state is: 
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Five differential equations are needed for the five unknowns: 
drop temperature, mass, velocity, number density and combustion gas 
velocity. The change in drop temperature is related to the heat received 
by the drop: 

DTt* _ Q* 

DJt* ‘VfCCfit 

The decrease of drop mass is equal to the vaporization rate: 

= - W* 

Pi"' 


(14) 


(15) 


The change of drop velocity is 

M = Cp z&Jl! 


2 


) (16) 


where Cj^ is the drag coefficient and is chosen to be 
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Re 


The number of droplets is conserved. 


so we have: 


(17) 


1 + v-n*|| = o 

The combustion gas continuity equation gives the relation for gas velocity 
and the mass addition: 



(18) 


Ifl + V p*^ = n* W* 


( 19 ) 


With proper initial and boundary conditions, this set of equations 
can be solved by numerical integration. 

In order to take a closer look at the effects of all the physical 
parameters, a dimensionless formulation is more valuable. The reference 
values of pressure and temperature are chosen to be the mean gas pressure 
and temperature, respectively - and the reference gas density is: 



, u 


Cf>^ , ju Oo > j ^ and 

The specific heat ratio is: 


CpiToo^^^ evaluated at 


if = Cp„/(Cpoo - rVm^ ) 


The speed of sound is: 

= j !r R-* T- /M| 


The reference length and velocity are chosen to be the chamber length 
and speed of sound, respectively and the reference time is defined as 


X 


Oo 


Lc^ 

Coo 


The initial drop radius is used to normalize the drop radius. 
Also, the reference quantities of the droplet > ^oo , and 

are evaluated at the initial temperature of the droplet. The initial 
drop number density is used to normalize the drop number density. 
Define the dimensionless quantities: 
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y= - 'T* le* 

T ’ 


I OO 


Cp^ 


OO 


f>< = 


(^O* 


_ZL 

A.OO 


For convenience, we separate Re into two parts: 


Re = Re^o Re 


where 


and 



Re OO — 2 Too Coc 

Re = ^ p 1 ^ ^ I y/ 


n = rr 


n 
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(24) 


Assume that: 


Sc = 
p; = 


Scoo ■= cotasi'atni" 
Pfoo ~ CO VI sto.iai’ 


We have Num and Nuh 

N um 


where 




= 2. 

0.6 4^ 

Re" 


= 2. 

+ o.G 'f" 
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Re^xa 
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(25) 

X 

X 



Reoo 

PrL 
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t y —I— * 

^ is related to by such relations: 

'f’v = exp (ol. -»• ) 


where ar . di -<• 

In dimensionless form; 

e^' 


where 


i’-.x = 


d. 
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<^o 


A. 

T. 


4i 

empirical < 

exp ( 

C^Z 
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> 
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Certain dimensional and dimensionless groups are defined as follows: 


( 27 ) 
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u*= in. 

= 4--K. x-i, 

W-r = "n. w 

The differential equations for drop temperature, drop mass, drop velocity, 
number density and gas velocity take the following forms; 


(33) 

(34) 
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(37) 

(38) 

(39) 


where c{ , (3 > S are dimensionless quantities which are 

fixed for a particular problem: 
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( 42 ) 



( 43 ) 


tts n«o -ji- 

1 

A numerical example of the normal heptane with gaseous oxygen case has 
been obtained for the steady mean flow; some of the physical quantities 
are given here: 

Too =- sooo.o 

L&o ■= 2-S.O 

Up 230.0 '^/sec- 

Tcso — 75 X io"‘^ 

Tc6o= 300*0 
U*„ = 30.0 ""/se^ 


Figures 34-38 show the variations of drop temperature, drop radius and 
mass, gas velocity and drop velocity, vaporization rate per drop and 
vaporization rate pi^r unit volume along the axial direction of the com- 
bustion chamber. Discussion of the results can be found in Ref. 28. 

By using the vaporization model just described, a quasi-steady linear 
analysis with transverse waves is developed. A small parameter £ is 
introduced here. This parameter is assumed to be the amplitude of pres- 
sure oscillation at a specific radial and angular position. 

We have 

P= l + fe p, 

(A ■= (A "V ^ CA , ^ U Q ^ ^ U -ft ^ 

'h' = € “T , ^ ^ , 


and so on. 


79 - 


It is assumed that ju, . p.D, , 4 , , 7\., . . Cpi, . 

can be neglected and also it is assumed that 
Assiime the isentropic relation for the gas 

i 

c 1 + e p, ) = C ^ ip' ^ 

C v+ t T. ■) = C fe -t>. ) 

Neglect the higher order terms, we obtain 

W ^ Vs 


With the assumption of pure transverse mode, (i.e., the 
velocity perturbation enters into the problem in higher order terms only. 
This makes the linear analysis possible. 


(44) 

(45) 


After some algebraic manipulation, the following perturbation quantities 
are obtained: 



(46) 


Mum I = 03 ^ Rc J (47) 

Mu4v \ ■=■ 0 <^ (48) 


Fi = V" , (49) 

to-. = _ Te , 

(Tc 


( 50 ) 



The perturbation of vaporization rate is 


W> = c, \£, -V -V 

where 
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2 - 
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[ Muw, + o.^ 4> 
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they are functions of the steady solution only. 

The perturbation of the vaporization rate per unit volume is 



Ugp 

Ug 


Other perturbation quantities are: 



Mum 



diL 

(Te-vclit 


^ C I - ^ ) 1,, 
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Only two partial differential equations remain to be considered here 
one for the perturbation drop mass and one for the perturbation drop 
temperature. They take the form of: 


^ + ue = A,T«, + A.'M, t Aj1p, 


= B,T«, i- > -V Si 
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(67) 


62 -=: -(SC. 

All > ^2 ’ ^3 ’ ^1 ’ ®2 ^3 functions of the steady- 

state solution only. 

A dimensionless frequency is introduced so that 


f = f 


( 68 ) 


(69) 


00 


(70) 


and the angular frequency is: 

CO == 2TTf 


( 71 ) 


The acoustic mode pressure oscillation has this form: 

-f, — P , A , G } .-c^ Cot 
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where P is a function of position. However, in the case of transverse 
oscillation and with the consideration of fixed radial and angular 
position, the pressure perturbation can be simplified to; 



P is absorbed in the 6 because it is a constant here. All other 
quantities will have the same transverse behavior since the velocity effect 
is of higher order. Therefore, only x and t dependence will be 
considered . 

We expect solution of these forms: 

T? 1 ■= "T<.s (73) 


<h\, = (W 

Thus, the space dependent part and time dependent part can be s epa rated 
and we obtain four ordinary differential equations for Tsis ’ > 

A ^ -V Qj) \ €c. 

(A:i 

AvT?c Ax^C - U3Ti2s (76) 

(Ax 

Cii! cWs _ + Bj -I- (-0 ‘Me oi) 

Ax 

Uii = Q>i Tcc -t — co^s (78) 

c\x 


This set of differential equations is ready to solve with proper initial 


conditions . 
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It is assxjmed that the droplet initial condition is not affected by 
the oscillation in the combustion chamber, therefore, the initial 
condition for the perturbation quantities S ’ 

can be taken as zero. 

So we have 


lesCo) = TzcCo) — "HsCo) ^ = o 


The vaporization rate perturbations are: 


W| = X Wc Cjtd CjO k (79) 

VMt.= Wts Co Jt -V Wt C Co-Cs CO “k ( 80 ) 

where 

Ws = C, Tij Cj. ITlj + 


\A/c ~ Uc “V- C-2. Wc. 


(82) 


W-rs = iAis 

IA/tc = ^ Uic 


(83) 


(84) 


Because numerical integration is required for the steady solution, the 
solution of the perturbation equations requires the result from the 
steady solution, therefore, numerical integration should be conducted 
simultaneously for both steady-state and perturbation equations. 

The idea of response factor is based on the Rayleigh criterion 
for acoustic amplification by heat or mass addition. Reinforcement of 
acoustic oscillation takes place when an excess amount of heat or mass 
is added while the pressure is above the mean value. The response factor 
is defined as the integral value of such mass addition over a given period 
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of time and in a finite volume normalized by the magnitude of the 
pressure perturbation. 


N ^ 


lA/r, i^i oiV^ 


^ Jri 

where ^ is the period of oscillation 


Jt = 


2.7c 

Co 


This response factor is positive when both vaporization rate per- 
turbation and pressure perturbation have the same signs and is negative 
when they have different signs. We can say that this response factor is 
an indication of the degree of driving or damping which the vaporization 
process gives a feedback to the system. 

At a specific radial and angular position the response is: 


/V = Jo Ia/ts 


A typical response curve is given in Figure 39 with the same 
physical condition as in that steady-state solution. 

At very low frequency, the response factor is nearly a positive 
constant. As the frequency increases, the response factor increases and 
reaches a peak value. Beyond that the response factor decreases and goes 
into the region of negative response. The critical frequency is defined 
when the response factor is zero. If the acoustic frequency of the com- 
bustor is in the range of positive response, particularly near the peak, 
the operation of such a combustor is most likely unstable. 

Several physical quantities which can affect the response factor 
are shown in Figures 40-45. 

Figure 40 shows the effect of combustor pressure. The higher the 
pressure, the larger positive response factor is obtained. Figure 41 
shows that the higher the combustor temperature the smaller the positive 
response factor. Figure 42 shows the effect of gas exit velocity. The 
effect on the peak is small but the critical frequency increases as the 
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gas exit velocity increases. Figure 43 shows that the higher the drop 
initial temperature, the lower the peak will become. This effect is 
important because the variation of drop initial temperature can be 
large. The effect of drop initial radius is shown in Figure 44 . 

Decrease of drop radius shifts the response curve to the right. Such 
effect is very important because the drop radius varies in quite a large 
range. Figure 45 shows the effect of drop initial velocity. The larger 
positive response corresponds to higher drop initial velocity. 

The result of linear analysis is similar to that of the nonlinear 
numerical analysis in Ref. 30 , as far as the response is concerned. The 
reason for these two different approaches giving approximately the same 
result is that the vaporization rate depends on the absolute velocity 
difference of the gas and the droplet. During a period of oscillation, 
the vaporization rate perturbation has two peak values and the effect of 
pressure gives nonlinear behavior of vaporization. However, such nonlinear 
effects will be canceled in the evaluation of response factor. Therefore, 
the velocity difference effect can be neglected as an approximation, that 
makes the linearization possible. 

In applying the vaporization mechanism as the combustion model, we 
must be aware of the limitations of the assumptions we have made. For 
instance, the uniform temperature over the whole droplet is quite question- 
able because the heat conductivity of the liquid droplet is not infinite 
and the intensity of the internal circulation of the droplet is not known. 

Let us examine the heat conduction equation for spherical symmetric 

case: 


il! = J_ ^ 

9r* 
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in dimensionless form: 


-feoo 
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(89) 



taken to be the period of oscillation here. The parameter 

\P0 has a dimension of time and is the characteristic 

time for heat diffusion. The uniform temperature assumption can be 
made only when 

<< \ 

We can see that this assumption is good only for small drop radius 
and low frequency oscillation. We can expect that the response factor 
we have obtained is not valid for large drop size and for high frequency 
oscillation. 

Finally, the response factor is not a criterion for combustor 
stability. A complete dynamic analysis of the whole combustion system 
is needed to acquire such stability information, further study in this 
direction will be conducted. 
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NOMENCLATURE 


A( , , Aj 

, B., B. 

c 

C.,C.,C3 

Cd 

Cp 

D 

C^l , J C\: 

f 

F 

H 

\ 

L 

TTl 

M 

Tl 

N 

Na\ 

(\lum 

tv 

P 

Pr 

t 

Q( 

Qv 


defined by Equation (64) -(66) 
defined by Equations (67) -69) 
speed of sound 

defined by Equations (53) -(55) 
drag coefficient 

specific heat at constant pressure 
diffusivity 

empirical constants in Equation (27) 
frequency 

defined by Equation (2) 
defined by Equation (3) 
t herma 1 c onduc t iv i t y 
length 
drop mass 
molecular weight 
drop number density 
response factor 

Nusselt number for heat transfer 

Nusselt number for mass transfer 

pressure 

vapor pressure 

defined by Equation (72) 

Prandtl number 
velocity 

defined by Equation (10) 
defined by Equation (6) 
radial coordinate 
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Y 

R 

Re 

Sc 

t 

T 

a 

vr 

\^/ 

Wt 

z 

c{ 

S' 

s 

e 

n 

0 

TV 



CO 


drop radius 
gas constants 
Reynolds number 
Schmidt number 
time 

period of oscillation 
temperature 

velocity component in axial direction 

velocity component in transverse direction 

vaporization rate per drop 

vaporization rate per unit volume 

axial distance from injector 

defined by Equation (9) 

correction factor 

defined by Equation (40) 

defined by Equation (41) 

specific heat ratio 

defined by Equation (42) 

amplitude of oscillation 

defined by Equation (43) 

angular coordinate 

latent heat 

viscosity 

density 

defined by Equation (25) 
defined by Equation (26) 


angular frequency 
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( L 

n 

( f 
( ). 


( )c 
( )s 

( )o 
1 ). 


reference value 
steady mean value 
dimensional value 
liquid or drop 
vapor 
gas 

coefficient of cosine component 
coefficient of sine component 
at injector face 
at nozzle entrance 
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D. NUMERICAL CALCULATION OF WAVE PHENOMENA IN CHAMBERS 

One of the Interesting questions about nonlinear gasdynamic phenom- 
ena in rocket chambers is whether or not shock waves may form in the 
transverse mode. Usually experimental evidence indicates continuous waves 
except in the longitudinal mode; however, in certain cases shocks have 
been observed. 

The problem of nonlinear transverse waves in a circular cylinder has 
been analyzed by Maslen and Moore who found solutions without shocks. 

There is, however, no indication that those solutions are unique. One way 
to determine if other solutions exist is to numerically integrate the 
equations of motion and see what the limiting behavior is after a long 
time. A variety of amplitudes and wave shapes could be used in the initial 
conditions and it would be seen if any of these resulted in the formation 
of shock waves. 

It is expected that shock waves would form in a thin circular annular 
chamber. Here the flow field is essentially one-dimensional (azimuthal) 
and centrifugal effects are negligible. One wonders how "thick” the 
annular chamber may be before shocks no longer form. An interesting param- 
eter, therefore, is the ratio of the inner radius to outer radius which 
varies between zero and unity. 

The problem remains as to how to write the difference equations to 
study transverse trav^slling waves. A popular technique for solving steady- 
state gasdynamic problems is to solve an initial-condition problem which 
yields the steady-state solution as its limiting behavior for large time. 
Note that for travelling waves a steady-state exists if one moves with the 
wave. This means that, if one were to observe the travelling wave from a 
coordinate system which rotates about the center of the cylindrical chamber 
at the same frequency as that of the oscillation, he would, after a long 
time, see a steady-state situation. 

It has been intended to employ this technique for studying the wave 
motion in the chamber. In order to insure rapid convergence and stability 
with this type of difference scheme, certain terms which have no physical 
meaning are added to the difference equations. In the steady-state limit 
the effect of these terms disappear and they only affect the transient 
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behavior. For this reason, the transient behavior is not accurately 
predicted by this method. 

In order to gain a better understanding of the advantages and 
disadvantages of this method, it was decided to compare it with the 
Lax-Wendroff scheme. The Lax method contains no artificial terms so 
that the transient behavior should be accurately predicted. Although 
the two-dimensional transverse case is of prime interest, it was decided 
to compare the two difference methods in the one-dimensional case for 
the initial attempt. Note that a travelling wave in a thin annulus is 
approximately a one-dimensional phenomenon, however. 

In these one-dimensional cases a steady-state solution could only 
be found if we observed from a frame of reference which was in transla- 
tional motion at the same velocity as the wave. Actually, in either the 
one or two dimensional cases, the continuous wave would distort until 
either a steady continuous waveform is obtained or a shock forms. If 
the shock forms no steady solution can really ever be achieved (in the 
absence of combustion) since dissipation occurs with a shock. After 
the shock is fully formed, its amplitude would decrease due to this dis- 
sipation. For weak shocks, however, this decay in amplitude is suffic- 
iently slow that a steady-state would be approximated. 

The equations of motion for the isentropic flow of a perfect 
gas written in a frame of reference moving at constant velocity k are 
as follows 



Here, the gas velocity u is measured with respect to the original frame 
of reference but the time derivatives are taken in the moving coordinate 

* This is only strictly valid if no shocks are present. Once shocks form, 
the isentropic assumption remains a good approximation as long as the 
shocks are weak. 
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system, p is pressure, yO is density, x is position, and t is 
time . 

The finite difference form of these equations are written as 

follows: 

nfl n ri.,n ^ \ P/ ^ 7 

/i “/i " ■ L -/ i-P “/ y i+i yOi-Pj 2/1. 

n f . r** X ^ n n x ^ n*l n*"l\ 7 t 

■/>i 1 *! - “i-i^ - r (“i+i - vi>J 
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^ t 
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1 r,, , r\/ n n . p, n-1 n-1. 7 At 

^ [a +r>(Pi+i - Pi.i) - / (Pi+1 - Pi.py 27^ 


where r^t and Ax are the step-sizes in the difference scheme, i and 
n are integers which denote the x and t positions, respectively, and 
r is a factor between zero and unity. 

Whenever P is zero, we have a simple difference scheme whereby 
the space derivative is represented by a centered difference. In other 
cases, that derivative is evaluated by averaging its value at two different 
times and p is a weighting factor in this averaging process. 

Also, the equations were written in conservation form in the orig 
inal coordinate system. Then, the Lax-Wendroff two-step difference scheme 
was employed. 

The same initial conditions were provided for both schemes. In 
particular, sinusoidal waveforms of linear theory were provided as initial 
conditions. In addition to waves travelling in one Cartesian direction, 
the case of travelling waves in an annulus was considered. In the latter 
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case, periodic conditions (in space) were applied. 

Figures 46 and 47 show pressure versus the azimuthal dimension 
in a thin annular chamber at fixed times. The first figure shows the 
results of the {" -method and the second figure shows the results of the 
Lax-Wendroff method. Both figures indicate a steepening of the waveform 
and a tendency towards shock formation. Qualitatively, the agreement 
between the two methods is good; there are, however, quantitative differences. 

It can be shown by means of Riemann invariants that, prior to 
shock formation, the minimum and maximum values of pressure will not change 
with time as the wave distorts. This is found with the P -method but not 
with the Lax method so that in this regard the -method is superior. 
Furthermore, it indicates that, for the present purposes, the Lax method 
may not be entirely reliable as a standard for comparison. 

No esgential differences were found between the processes of 
shock formation for waves travelling in one Cartesian direction and waves 
travelling in a thin annulus. 

It was concluded that the f'' -method provides a good prediction 
of the wave phenomenon even for the transient behavior. The difference 
equations were not written in a form which could account for shock waves so 
that a steady-state was never actually achieved. Since the ultimate purpose 
of this endeavor is to determine whether shocks form in the two-dimensional 
transverse case, it was decided not to spend any time to include the capability 
of handling discontinuities in the computer program until their formation 
became evident in the two-dimensional case. (Of course, in that event, a 
technique such as the pseudo-viscosity method may be employed.) For this 
reason, the one-dimensional calculations ended after considerable distortion 
of the wave but before discontinuities formed. 

This endeavor is continuing with the programming of the two-dimen- 
sional case in annular and "full" chambers with both difference schemes. 
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